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ABSTRACT 

We examine the K shell emission lines produced by isothermal and simple multi- 
phase models of the hot gas in elliptical galaxies and galaxy clusters to determine the 
most effective means for constraining the width of the differential emission measure 
(£(T)) in these systems which we characterize by a dimensionless parameter, er^. Com- 
parison of line ratios of two-temperature (cr^ <C 1) and cooling flow (cj^ ~ 1) models 
is presented in detail. We find that a two-temperature model can approximate very 
accurately a cooling flow spectrum over 0.5-10 keV. 

We have re-analyzed the ASCA spectra of three of the brightest galaxy clusters 
to assess the evidence for multiphase gas in their cores: M87 (Virgo), the Centaurus 
cluster, and the Perseus cluster. Kct emission line blends of Si, S, Ar, Ca, and Fe are 
detected in each system as is significant Fe K(3 emission. The Fe K/3/Ka ratios are 
consistent with optically thin plasma models and do not suggest resonance scattering 
in these systems. Consideration of both the ratios of H- like to He-like Ka lines and the 
local continuum temperatures clearly rules out isothermal gas in each case. To obtain 
more detailed constraints we fitted plasma models over 1.6-9 keV where the emission 
is dominated by these K shell lines and by continuum. In each case the ASCA spectra 
cannot determine whether the gas emits at only two temperatures or over a continuous 
range of temperatures as expected for a cooling flow. The metal abundances are near 
solar for all of the multiphase models. We discuss the implications of these results and 
examine the prospects for determining the temperature structure in these systems 
with upcoming X-ray missions. 
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1 INTRODUCTION 

There is mounting evidence that the X-ray emission in el- 
liptical galaxies and galaxy clusters originates from multi- 
phase plasma. The first important evidence for multiphase 



gas in e llipticals and clusters came from analysis of emission 
lines in M87 and the Perseus cluster. Using ratios of sev- 
eral blends of Fe L shell lines obtained from the Focal Plane 
Crystal Spectro meter (FPCS) on the Einstein Observatory, 
Canizares et al. ([L98§) found that the hot gas in M87 emits 
over a large range of temperatures: from ~ 2 — 4 x 10 6 K to 
~ 3 x 10 7 K. The strength of the Fe XVII line measured from 
FPCS data of the Perseus cluster also suggests a multiphase 
plasma (Canizares, Markert, & Donahue 1988). 



ASCA (Tanaka et al. 1994) satellites have provided new ev- 
idence for multiphase hot gas. Broad-band spectral analysis 
of ASCA data of clusters (Fukazawa et al. 1994; Fabian et al. 
1994a), the brightest ellipticals (Buote & Fabian 1998; Buote 
1999a) , and of the brightest poor galaxy groups (Buote 



1999b) indicates that the hot gas in the central regions of 



these systems consists of at least two temperature compo- 
nents. The equivalent widths and ratios of Ka line blends 
of Si and S alone strongly favor multitempe rature model s of 
the ASCA data of the brightest ellipticals (Buote 1999a). 



Recently, data from the ROSAT (Trimmer 1983) and 



Complementary evidence for multiphase gas in clusters 
i s provided by the analysis of imaging data from ROSAT 
(Trimmer 1983) of cooling flow clusters which infers mass 
deposition rates which vary with radius approximately as 
M ex r (e.g. Peres et al. 1998) . This type of distributed mass 
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deposition is the signature of multiphase cooling flows (for 
a review see Fabian 1994). Multiphase cooling flows are also 
consistent with the positive temperature gradients observed 
in the centers of the brightest elliptical galaxies with ROSAT 
(e.g. Forman et al. 1993; David et al. 1994; Trinchieri et al. 
1994; Rangarajan et al. 1995; Irwin & Sarazin 1996; Jones 
et al. 1997; Buote 1999a). 

It is essential that the temperature structure be ac- 
counted for in order to deduce correctly the properties of 
the hot gas in these systems. The metal abundances are es- 
pecially sensitive to the assumed temperature structure. For 
example, analyses of the X-ray spectra of elliptical galaxies 
where the hot gas is assumed to be isothermal always ob- 
tain very sub-solar metallicities, Z ^ O.SZq, whereas when 
multi-temperature models are considered the metallicity of 
ellipticals becomes consistent with solar (see Buote & Fabian 
1998 and Buote 199 9a). A simila r situation prevails for poor 
groups of galaxies ( Buote 1999b| ). 

Although this "Fe bias" is unimportant for clusters 
of galaxies with temperatures in excess of ~ 3 keV, if an 
isothermal gas is assumed for a cooling-flow cluster (or el- 
liptical or group) then the Si/Fe and S/Fe ratios can be sig- 
nificantly over-estimated (see the appendix in Buote 1999b) . 
Hence, in order to obtain reliable measurements of the heavy 
element abundances crucial to interpreting models of the for- 
mation, evolution, and enrichment of the hot gas in ellipti- 
cals, groups, and clusters (e.g. David, Forman, & Jones 1991; 
Loewenstein & Mathews 1991; Ciotti et al. 1991; Loewen- 
stein & Mushotzky 1996; Renzini 1997; Brighenti & Math- 
ews 1999), it is necessary to know the true temperature 
structure of the hot gas. 

Knowledge of the temperature structure is also re- 
quired to understand the contribution of discrete sources 
to the hard X-ray emission in elliptical galaxies. For exam- 
ple, if the hot gas is assumed to be isothermal then the 
hard spectral component measured by ASCA for the ellip- 
tical galaxy NGC 3923 is inconsistent with the expectation 
that discrete sources follow the same profile as the optical 
light (Buote & Canizares 1998; see also Pellegrini 1999). 
Finally, when multiphase cooling flows are accounted for 
in X-ray observations of clusters the scatter in the X-ray 
luminosity-temperature relationship is significantly reduced 
(e.g. Fabian et al. 1994b). 

The Fe L shell emission lines from ~ 0.7 - 1.4 keV gen- 
erally provide the strongest constraints on the temperature 
structure of ellipticals, groups, and clusters. There has been 
considerable progress in improving the accuracy of the the- 
oretically computed Fe L shell transitions (e.g Liedahl et al. 
1995), but they are still among the most uncertain transi- 
tions implemented in the currently available plasma codes 
(e.g. Savin et al. 1999). There is reason for some confidence 
in the codes since recent empirical studies have concluded 
that the errors in the Fe L transitions are not overly serious 
(Hwang et al. 1997; Buote & Fabian 1998). (However, see 
the recent papers by Savin 1999 and Savin et al. 1999 for 
discussions of the inaccuracies in the dielectronic recombi- 
nation rates used by the popular plasma codes which affect 
all types of emission lines.) 

Nevertheless, it is important to use lines from the more- 
accurate K shell transitions to search for evidence of multi- 
phase gas to complement the information from Fe L lines. 
To obtain the most accurate constraints on the emission 



measure of a multiphase plasma it is necessary to obtain 
independent measurements of the temperature(s) over the 
entire energy range of significant emission. It is a happy co- 
incidence that the K shell lines of many abundant elements 
appear throughout the X-ray wave band (0.1-10 keV). 

Line ratios have important advantages over broad-band 
spectral fitting for probing the temperature structures of el- 
lipticals and clusters. Firstly, in an isothermal plasma using 
ratios of lines of the same element avoids uncertainties in the 
relative abundances of the elements; there is a small depen- 
dence on metallicity in a multiphase cooling flow (section 
2.2.3). Secondly, ratios of lines nearby in energy are essen- 



tially independent of Galactic or any intrinsic absorption. 
Thirdly, since the strength of a line can be computed using 
a local continuum measurement, line ratios are insensitive to 
any poorly constrained power-law or bremsstrahlung spec- 
tral components. 

For a multiphase coronal plasma Q the luminosity of an 
emission line arising from a transition from an upper state 
u to a lower state I of an element X in ionization state r 
may be expressed as, 



L ul {X r ) = / A ul (X r ,T)£(T)dT, 



(1) 



where A u i(X r , T) is the plasma emissivity of the line at tem- 
perature T which incorporates the abundance of element X, 
the relative fraction of ionization state r, and other quan- 
tities depending only on the atomic physics. (The integral 
extends over the emitting volume V .) The quantity, 

£(T) = n*(T)^, (2) 

is called the "differential emission measure" and £(T)dT the 
"emission measure" at temperature T. The emission mea- 
sure weights A as a function of temperature in equation ([!]) 
and thus effectively characterizes the temperature structure 
of the plasma. Notice that different phases T in general can 

have different electron number densiti es n e . 

It has been known for some time dCraig fe Brown 1976; ) 
that precise constraints on £(T) from inversion of equation 
(^) are difficult to obtain in practice. (For a recent discussion 
of this topic see, e.g., Judge, Hubeny, & Brown 1997.) That 
is, for a given line typically A u i(X r ,T) is significant over 
temperature ranges of factors of ~ 2-5, which means that 
the plasma emissivity acts as a smoothing function which 
obscures the detailed shape of £(T). Moreover, small un- 
certainties in the data are amplified by the inversion. Given 
these difficulties, it may only be possible to obtain relatively 
crude constraints on the shape of £(T)dT for ellipticals and 
clusters from inversion of equation (ill); e.g., width, skewness, 
etc J} 

The primary theoretical goal of this paper is to iden- 
tify the best X-ray emission line ratios for determining the 
gross features of £(T) of ellipticals and clusters with the spe- 
cific objective of distinguishing between single-temperature, 
two-temperature, and cooling-flow spectra. (Previous dis- 



* The assumptio ns of the standar d coro nal model are discussed 
in, e.g., Sarazin (L98£) and Mewc (1991). 

T See, e.g., Kaastra et al. (1996) and Schmitt et al. (1996) for 



recent constraints on §(T) of stellar coronae. 
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cussions of important temperature sensitive lines for clus- 
ters in the context of a diffe rent set of models m ay be found 
in Bahcall & Sar azin ( 1977), Sarazin & Bahcall (1977), and 
Canizares et al. (1982). Wise & Sarazin (1993) discuss the 
radial variation of the shapes of emission lines as diagnos- 
tics for cooling flow models.) To facilitate quantitative com- 
parison of the multiphase structure of different ellipticals 
and clusters we introduce a quantity, erg, the "multiphase 
strength" . We focus on K shell emission lines since they are 
generally the strongest lines of ellipticals and clusters and 
the atomic physics is relatively simple and so the available 
plasma codes should be most accurate for these transitions. 
Because the Fe L shell lines occurring at energies ~ 0.7— 1.4 
keV are strong for ellipticals and clusters and are powerful as 
temperature diagnostics, we also examine them even though 
at the present time the accuracy of the plasma codes is less 
certain for these transitions (see above). The line ratios allow 
simple yet effective discrimination between popular spectral 
models of ellipticals and clusters and serve as an important 
guide to interpreting results of £(T) obtained from inversion 
of equation (hi). 

The paper is organized as follows. In Section B we ana- 
lyze theoretical line ratios in single-phase (Section |2.l[ ) and 
multiphase (Section 2.1) mode ls. The multiphase strength 
is introduced in Section 



2.2.1 



and it is used to elucidate 
the tem perat ure structure of simple multiphase mod els in 
Section 2.2.2 and of cooling flow models in Section 2. 2. J . 
We discuss in detail how to distinguis h a two-temperature 
plasma from a cooling flow in Section 2.3. In Section H we 



compute and analyze line ratios of ASCA data of M87, the 
Centaurus cluster, and the Perseus cluster. Simulated obser- 
vations with Chandra, XMM, and ASTRO-E are considered 
in Section ^. Finally, we comment on recent evidence for 
resonance scattering in 5.1 and discuss cosmological appli- 
cation of the multiphase strengths of clusters in Section [5.2| 
In Section ^ we present our conclusions. 



2 THEORETICAL LINE RATIOS OF 
ELLIPTICALS AND CLUSTERS 

We examine the theoretical emission lines predicted by coro- 
nal plasma models of ellipticals, groups, and clusters. In or- 
der to maintain a connection to feasible observations within 
the next 10 years, for the purposes of identifying lines we 
consider spectra folded through an energy resolution of 2 
eV which is the targeted performance of the Constellation- 
X mission To model the emission of a coronal plasma we 
use the MEKAL code which is a modification of the original 
MEKA code (Mewe, Gronenschild, & van den Oord 1985; 
Kaastra & Mewe 19 93) where the Fe L shell transitions have 
been re-calculated (Liedahl et. al 1995). We fix the relative 
abu ndances of the elements to the (photospheric) solar val- 
ues (Anders & Grevesse 1989) and, as a result, we focus on 



ratios of lines of the same element so that our results are 
not sensitive to this assumption. Finally, we do not include 
photo-electric absorption (Galactic or intrinsic) in this dis- 



t http:/ /constellation. gsfc.nasa.gov/ 
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Table 1. K Shell Transitions for H-like and He-like Ions 

Isoclcctronic Transition 
Sequence Symbol 

H 1 ls 2 Si 

2 
3 



Is 2 Si 
Is 2 S? 



Transition 

- 2p 2 Pi 

- 3p 2 Pj 



He 



Is 2 Si 

Is 2 !&} 

Is 2 % 
Is 2 

Is 2 % 

Is 2 l So 

Is 2 % 



4p Ti 3 

„ 3-7 

5p 2 Pl 3 

2 ' 2 

ls5p x Pi 
ls4p x Pi 
ls3p 1 Pi 
ls2p J Pi 
ls2p 3 P 2 , 
ls2s 3 Si ' 



Transitions are given in the standard notation of the electron 
configuration (orbitals) and the level with the lower state listed 
first. HI and He4-6 are all K« transitions while H2 and He3 are 
K/3 transitions. 

cussion^] but include its effects when interpreting line ratios 
from real data (as in Section []). 

The emission lines which are perhaps the most useful 
for our study of the temperature structure of ellipticals and 
clusters pj| are the K shell transitions involving ions of the 
H-like and He-like isoelectronic sequences listed in Table ^ 
using the standard notation (e.g. Table 3 of Mewe et al. 
1985). Also very useful are the satellites to the resonance line 
of the preceding ion in the ionization sequence. Of particular 
interest are the satellite lines to the He4 resonance from 
dielectronic recombination to autoionizing states of the Li- 
like ion. The strongest of these satellite lines typically come 
from the multiplet ls 2 2p 2 P - ls2p 2 2 D. Finally, the L shell 
transitions of Fe between ~ 0.7 - 1.4 keV (especially the 2s- 
3p lines) are strong and very temperature sensitive, though 
most are unresolved at 2 eV resolution and thus blends of 
lines are emphasized below. 

We identified significant K and Fe L shell lines between 
0.4 - 10 keV from visual inspection of isothermal MEKAL 
models. These models were viewed at 2 eV resolution and 
had temperatures ranging from 0.1 - 10 keVjJJj We restricted 
our attention to this energy range because below ~ 0.4 keV 
Galactic photo-electric absorption is prohibitive, and ener- 
gies above ~ 10 keV will not be effectively probed by high 
energy resolution X-ray satellites expected to operate in the 
near future. 

For lines that are not isolated at the 2 eV resolution, 
our chosen energy boundaries in many cases enclose other 
emission lines from the same element or others. The result 
is that the line ratios we discuss below often behave differ- 
ently over some temperature ranges than expected if there 
were no blending of lines. However, we enforce that such pe- 
culiarities are not the result of uncertainties in defining the 
continuum because the line luminosities are computed from 
the models using a finer energy resolution which allows the 



S We di scu ss the possible complication of resonance scattering in 
Section [a.l[ 

» A summary of standard tem perat ure diagnostics of coronal 
plasmas may be found in Mewe (1991). 



II Throughout the paper we quote temperatures in units of keV 
(k B T). 
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Table 2. Energy Ranges of Satellite Lines 



Ion(s) 
Si XII 
S XIV 
Ar XVI 

Ca XVII/XVIII 
Fe XX/XXI 
Fe XVIII/XIX 
Ni XXV 



Energy Range 
(keV) 
1.847-1.852 
2.436-2.445 
3.130-3.137 
3.863-3.875 
6.610-6.630 
6.520-6.610 
7.725-7.747 



Isoelectronic 
Sequence 
Li 
Li 
Li 

Be/Li blend 
N/C blend (satl) 
F/O blend (sat2) 
Li 



Energy ranges used to define the dielectronic recombination satel- 
lite lines. 



continuum to be accurately subtracted. Within the energy 
boundaries used to define a line (or blend) we compute the 
total emission above the continuum on an energy grid of 
15000 logarithmically spaced energy bins over 0.4 - 10 keV; 
i.e. ~ 0.2 eV bins around 1 keV and ~ 1.5 eV bins around 
7 keV. 

Our procedure to define the continuum for the lines 
of interest begins by identifying from visual examination of 
the spectrum a nearby bin at lower energy having no ap- 
parent lines. This choice is made after examining spectra of 
several temperatures. (Note that we avoid energy bins hav- 
ing strong recombination edges.) Let us denote the contin- 
uum emission at this energy bin, j c (Ei), where Ei is the 
energy of this bin i. Proceeding to the next bin (i + 1) 
with energy -Ei+i, if there is no line emission in this bin, 
then the total emission j(Ei+i) is all continuum such that 
j(Ei+i) < jc(Ei) since (at a fixed temperature) the (free- 
free) continuum is a monotonically decreasing function of 
energy. In this case, j c (Ei+i) = j(Ei+i) becomes our new 
definition of the continuum. If there is line emission in bin 
i+1 so that j(E i+1 ) > jc(E»), then we set j c (E i+1 ) = j c (Ei). 
This procedure is followed up to and through the lower 
boundary of the emission line of interest until the upper 
boundary of the line region is reached. The continuum values 
(jc(Ei),j c (Ei + i), . . .) are subtracted from the corresponding 
bins of the line region. 



2.1 Isothermal case 

We begin by examining isothermal plasma models in detail 
since the temperature dependences of line strengths and line 
ratios are most easily understood for these models. We plot 
spectral regions in Figure [l] corresponding to the strongest K 
shell transitions of Si, S, Ar, and Fe. In Figure^ we plot the 
luminosities of the K shell emission lines of Si XII-XIV and 
Fe XXIII-XXVI obtained from the isothermal models; these 
ions have among the strongest K shell lines which are also 
isolated from strong lines of other elements. Since the Li-like 
satellites of Fe XXIV are unresolved at 2 eV resolution, the 
Fe satellites plotted in Figure ^ result from ions in the next 
few isoelectronic sequences. In Table ^| we list the energy 
ranges used to define these satellite blends for Fe as well as 
the Li-like satellites of the other elements. 

In Figure ^ we plot as a function of temperature the ra- 
tios of the K shell emission lines of the abundant elements. 
(Ratios involving lines that are very contaminated are not 
plotted.) The ratios of lines emitted from subsequent ion- 
ization stages of the same element (Hl/He4, Satellite/He4) 



are overall more useful as temperature diagnostics than are 
the ratios involving lines of the same ion (H2/H1, He3/He4). 
The former, however, rely on the assumption of ionization 
equilibrium which can be checked using the latter ratios. 
Ionization equilibrium is expected to be an accurate approx- 
imation for the hot gas in elli pticals, groups, and c lusters, 
even those with cooling flows ( [Canizarcs et al. 1988j ) . 

The Hl/He4 ratio is probably the most useful because 
the HI and He4 transitions give rise to the strongest lines 
and the ratio ~ 1 where the lines have maximum luminosity; 
see the examples of Si and Fe in Figure]^. The mostly smooth 
curves of the Hl/He4 ratios as a function of temperature in 
Figure ^| demonstrate that the HI and He4 lines of many 
of the abundant elements can be measured accurately at 2 
eV resolution. The largest distortions are evident for the Mg 
ratio which result from contamination in the measurement of 
the Mg XI He4 transition by Fe L emission. Other significant 
contamination appears in the computation of the Fe XXV 
He4 transition which includes unresolved Li satellites at 2 
eV resolution. 

Since the Hl/He4 ratios for the different elements take 
a value of 1 approximately uniformly distributed over the 
logarithmic range of temperatures ~ 0.5 - 10 keV, accu- 
rate estimation of any temperature over this range can be 
obtained by Hl/He4 ratios of more than one element. For 
temperatures below ~ 0.5 keV the only useful Hl/He4 ratio 
is that of oxygen. 

The next most useful diagnostic is the Satellite/He4 ra- 
tio, where "Satellite" typically refers to the Li-like satellites; 
see Table ^ for the definitions of the satellite lines used. The 
principal difficulty here is measuring the satellite lines be- 
cause they are usually much weaker than the He4 line and 
are difficult to resolve from the (forbidden) He5 line, espe- 
cially for elements with atomic number 20. The shapes of 
the temperature profiles of the Satellite/He4 ratios of Si, S, 
Ca, and Ni are similar in Figure |§| The ratio of Fe K lines 
has a steeper profile because satellites from isoelectronic se- 
quences above Li are used; i.e. satsl in Table ^[ The distorted 
temperature profile of the Ar ratio for T ^ 1 keV arises from 
unresolved Li satellites being included in the He4 line. 

Over T ~ 0.5 - 5 keV the Si Satellite/He4 ratio changes 
by a factor ~ 20 compared to ~ 1000 for the Si Hl/He4 
ratio. The steepest Satellite/He4 profile shown is for Fe K. 
Even if we do not consider T < 2 keV where the Fe K lines 
are weak, then the Satellite/He4 ratio changes by a factor 
of ~ 50 from T ~ 2 - 8 keV. This variation is about as large 
as seen for the Fe K Hl/He4 ratio. 

The H2/H1 (i.e. Ly /3/Ly a) ratios are much less 
temperature sensitive than are the ratios described previ- 
ously. The temperature dependence of this ratio is simply, 
H2/H1 oc exp[-(£2 - El)/k B T], where k B is Boltzmann's 
constant and E2 and El are the energies of the H2 and 
HI lines respectively. For Si, E2 — El — 0.37 keV, and as is 
clear from Figure || H2/H1 only varies significantly for tem- 
peratures within approximately a factor of 2 of 0.37 keV. 
At these temperatures the hydrogenic lines of Si are also 
very weak. This behavior holds qualitatively for the other 
elements shown. (Note that the plotted S ratio is actually 
H3/H1 because the H2 transition is blended with the Ar 
XVIII He6 intercombination line.) Typically, for tempera- 
tures where the lines are at their strongest, the H2/H1 ratio 
varies by at most a factor of 2; e.g. Si from 0.5-5 keV. 
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Figure 1. Each panel shows a portion of the X-ray spectrum where the K shell emission lines are prominent for an important abundant 
element: Si XII-XIV, S XIV-XVI, Ar XVII-XVIII, and (mainly) Fe XXV-XXVI. The spectra are computed for isothermal MEKAL 
models with Z = 0.5Zq and are plotted at 2 eV energy resolution; the temperatures are 1 keV for the Si and S panels, 2 keV for Ar, 
and 4 keV for Fe. The transitions mentioned in the text are labeled in each figure according to the notation in Table Be the labels include 
the name of the element only if it is not the element labeled on the panel. The luminosity is calculated for a total emission measure 
rtg V = 10 65 cm -3 appropriate for M87 (e.g Canizares et al. 1982) which is approximately midway between the luminosities of ellipticals 
and clusters. 




Similar to H2/H1, the ratio He3/He4 (i.e. K/3/Ka) in- 
volves lines of the same isoelectronic sequence and has a 
fairly weak temperature dependence for temperatures where 
the lines are of significant strength. Taking S for an example 
(Figure for temperatures ranging from 0.5 - 5 keV the 
He3/He4 ratio varies by only ~ 5 - 10 per cent. For T £ 0.5 



keV, the ratio is very sensitive to temperature, but the lines 
are also very weak: at T = 0.25 keV we have log 10 L x = 36.1 
and 36.7 for He3 and He4 of S assuming a total emission 
measure of 10 65 cm -3 (as in Figure ^). Probably Fe XXV 
is the most useful ion for this ratio because at ~ 2 keV 
where He3/He4 begins to change noticeably, the emission 
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lines still have significant emissivity; e.g., at T = 2 keV we 
have log 10 L x = 39.0 and 39.6 for He3 and He4 of Fe (see 
Figure 0). Hence, for Fe the He3/He4 ratios varies by about 
a factor of 10 over interesting temperatures. (Note that the 
Si ratio turns over below 0.2 keV because of weak Al XI 
contamination in the He4 line of Si XIII.) 

For any individual element, the ratios described above 
involve lines which differ in energy by < 1 keV. That is, each 
ratio gives a measure of the temperature in a small energy 
region of the spectrum. Comparing these ratios for different 
elements over the entire X-ray spectrum allows one to exam- 
ine the consistency of the temperatures measured by each 
ratio. Another method to compare lines that are substan- 
tially different in energy is to compare ratios of lines of the 
same isoelectronic sequences of different elements (Figure^). 
Although these ratios depend on the relative abundances of 
the elements, they can provide an interesting consistency 
check of the other ratios. (In fact, these ratios are used to 
measure abundances once £(T) is known) Typically, these 
ratios vary by a factor of 100 over the temperature ranges 
where the lines emit significantly (Figure ft). (These ratios 



have a stronger dependence on Nu since the lines are not 
close together in energy, but for energies above ~ 2 keV the 
dependence should be negligible.) 

Finally, we consider the Fe L shell emission lines. Even 
at 2 eV resolution most of the Fe L lines between 0.7-1.4 
keV are not resolved. From visual inspection of the MEKAL 
model spectra we identified blends of nearby lines which var- 
ied similarly with temperature and were least contaminated 
by lines of other elements. In Table ^ we list energy ranges 
for five blends (or, more accurately, groups) of lines and 
plot their strengths and ratios as a function of temperature 
in Figure ^. These line groups of ~ 10 - 20 eV width are 
shown because they have the most regular behavior over 
the temperature range investigated. (The lack of blends be- 
tween 1.1-1.3 keV is due to the difficulty of finding groups 
uncontaminated by emission lines of other elements.) It is 
clear from comparison to the K shell strengths and ratios 
discussed previously that the Fe L lines are very strong and 
very temperature sensitive over energies ~ 0.2 - 5 keV; e.g. 
the 2:1 ratio varies by a factor of ~ 1000 over 0.2 - 1.5 keV. 
The superior performance of the Fe L line ratios highlights 
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Figure 4. Ratios of K shell emission lines from the same iso- 
electronic sequences of different elements for isothermal MEKAL 
models with Z = 0.5Zq. 



Table 3. Energy Ranges of Fe L Blends 



Blend 
1 

2 
3 
4 
5 



Energy Range 
(keV) 
0.722-0.741 
0.867-0.878 
0.940-0.975 
1.375-1.387 
1.435-1.445 



Dominant Ion(s) 

Fe XVII 
Fe XVII & XVIII 
Fe XVII & XVIII 
Fe XXII 
Fe XXIII 



the importance of improving the accuracy of the plasma 
codes for these transitions. It may be that future observa- 
tions will shed light on the accuracy of the codes in the Fe 
L region by comparing to the temperatures inferred from K 
shell transitions. 



2.2 Multiphase cases 

We now examine the ability of the line ratios to quantify 
the general multiphase structure of the hot gas in ellip- 
ticals, groups, and clusters. The temperature structure is 
fully characterized by the distribution of emission measure 
as a function of temperature (see Introduction). Of particu- 
lar relevance to ellipticals, groups, and clusters are emission 
measures of a single temperature component, a cooling flow 
component (see below), or a mixture of both. The single 
temperature component is often modeled as a 5-function dif- 
ferential emission measure while a multiphase cooling flow 
emits over a wide range of temperatures. Hence, to facilitate 
our discussion of the temperature structures of ellipticals 
and clusters, we wish to define a dimensionless parameter, 
o £ , which isolates the information contained in £(T) indi- 
cating the multiphase "width" or "strength" . (We discuss a 
possible cosmological application of in Section [5.2j,) 



2.2.1 Multiphase strength 

We first consider the case where £(T) does not contain any 
5-function components. This is the physical situation since 
in a real elliptical or cluster the hot gas will always have 
some temperature fluctuations about the mean as a result 
of incomplete relaxation. 

Over some temperature range (T m i n , r max ) it is neces- 
sary to define the "maximal multiphase strength" which we 
take to be £(T) = constant; i.e. for T within this range the 
emission measure, £(T)dT, is constant for fixed intervals dT. 
For constant £ the total emission measure equals £AT, where 
AT = T max — T m in. If £ (T) is centrally peaked (i.e. approxi- 
mately a single phase) then £ max AT will greatly exceed the 
total integrated emission measure. Thus, a sensible defini- 
tion for the "multiphase strength" of an emission measure 
distribution is, 



jT^ aT)dT ^ ^ 

2(T> 



^max AT 



2<T)£„ 



£(T) dT, 



(3) 



(4) 



where £ max is the maximum value of £, (T) is the emission 
measure weighted temperature, 



(T) 



£™~TgT)dT 
jT m ^ aT)dT • 



(•») 



and (T m in, Tmax) are any two temperatures which bracket 
the region in temperature space where £ 7^ 0. The fac- 
tor AT/2(T) in equation ^ takes into account the size 
of the temperature interval; i.e. for uniform £ we have 
(j c = AT/2(T) « 1 for a narrow interval T max » T min 
(i.e. approximately single phase) but er^ = AT/2(T) ~ 1 for 
Tmax 2> Tmin- Thus, <7j is a measure of the fractional width 

ofe 

To illustrate the meaning of cr^ let us for the moment 
consider £ to be a Gaussian with standard deviation a and 
mean temperature (T). In this case there is an analytic ex- 
pression for the multiphase strength, 



2 (T> 



(Gaussian), 



(6) 



which implicitly assumes (T m i n ,T max ) = (— 00, +00) but is 
accurate provided AT ^> a. As would be expected, oc a, 
indicating that measures the width of the emission mea- 
sure distribution, but 05 applies equally well to ellipticals 
and clusters because the mean temperature, (T), is divided 
out. 

The definition of multiphase strength just described is 
applicable for all "well behaved" £(T). However, it is often 
convenient to fit models consisting of individual 5-function 
temperature components to X-ray data of ellipticals and 
clusters. This widespread practice is well justified for nearly 
isothermal systems because the true finite width of £(T) 
cannot be distinguished from a 5-function due to the limita- 
tions in the energy resolution of current X-ray data. Hence, 
since (^-function temperature components (even though they 
are not quite physical) are often fitted to data of ellipticals, 
groups, and clusters, we consider how to modify the defi- 
nition of <j£ so that it applies for emission measure distri- 
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butions having only 5-function temperature components or 
having both a 5-function and "well behaved" temperature 
components. 

For £(T) consisting of a single ^-function temperature 
component we have <j£ = because £ ma x = oo. This in- 
dicates correctly that the width of the emission measure 
distribution is zero. However, if we consider £(T) to be 
composed of a single 5-function component plus some "well 
behaved" component such as a Gaussian, then <tj is still 
zero regardless of the shape of the "well behaved" com- 
ponent. When considering emission measure distributions 
composed of both a 5-function temperature component and 
other non-singular components, it is reasonable to modify 
0£ according to the fractional contribution of the "well be- 
haved" component to the total integrated emission measure. 
That is, if EM ns and EMj are the total emission measures 
of the non-singular and 5-function components respectively, 
then the total multiphase strength is constructed by re- 
placing of the non-singular component by /erj, where 
/ = EM„ S /(EM ns +EM 4 ). This form of the multiphase 
strength is especially useful when considering models con- 
sisting o f a mi xture of a coo ling flow and ambient gas (as in 



Sections 2.2.3, 2.3, and 5.2) 



The definition of cannot be so easily modified if f(T) 
consists of multiple (^-functions. A differential emission mea- 
sure consisting of a finite sum of N 5-function temperatures, 
f (T) = y^- —1 £,iS(T — Ti), always has = independent of 
the £j and N. Clearly should depend on TV and on the 
temperature difference between each (5-function. To simplify 
let us restrict ourselves to the case where the (^-functions 
are evenly spaced in temperature. In the limit of large N we 
can generalize to the discrete case cr^ by simply replacing 
£ max AT in equation Q with iV£ max . This is appropriate for 
N 3> N c , where N c is some critical number of S- functions 
above which a discrete number of phases is a good approx- 
imation to a continuous £(T). In the other limit we require 
that cr^ — + as TV — > 1. Considering these limits we define 
o? as follows, 



A^max 



AT 

W) 



(1 



-(JV-1)/JV C 



(7) 



The first two terms are just the straightforward generaliza- 
tion of <7£ for large TV while the last term is an arbitrary 
means to smoothly connect the limiting cases. A reasonable 
choice for N c is ~ 20. 



2.2.2 Simple examples 

We plot the Hl/He4 ratio of Si as a function of tempera- 
ture in Figure [] for MEKAL models consisting of two tem- 
peratures with equal emission measures. For clarity we ex- 
press a™ as the product of two factors: (1) the terms which 
are the generalization of for large N, and (2) the fac- 
tor involving the number of temperature components which 
is 0.05 for N = 2 and jV e = 20. Even for relatively small 
<j£ = 0.25 x 0.05 there are considerable departures from 
the isothermal case. Of particular interest is that there are 
values of (T) where very different a™ have very similar line 



ratios. Just above 1 keV the ratios for models with o> 







and 0.5 x 0.05 differ by only ~ 15 per cent. Just below 
(T) = 3 keV the curves of af = and 0.9 x 0.05 intersect! 

The key implication of these curves is that one line ra- 
tio does not uniquely determine a™ . The situation is even 
more complicated than shown in Figure ^| if the two tem- 
perature components are allowed to have different emission 
measures. In that case models having the same (T) and 
can also have different line ratios. These degeneracies are 
related to the ambiguity in determining f (T) from inversion 
of equation (|^) (see Introduction). To narrow the region of 
a™ space, one must add constraints from similar ratios of 
other elements and use as many temperature sensitive ratios 
as possible (e.g. satellite/He4). 

The situation is more tractable for a simple non-singular 
differential emission measure such as a constant or a Gaus- 
sian. In Figure ^| we show the Si Hl/He4 ratio versus tem- 
perature for the case £ = constant. The ratio profile changes 
slowly away from the isothermal curve and flattens for larger 
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Figure 6. Hl/He4 ratio for Si versus emission measure weighted temperature for two-temperature models, f(T) <x 5(T — T\) + <5(T — 
T2) (left), and models with a constant differential emission measure, £ = constant, over some temperature range (right). Each curve 
corresponds to a different multiphase strength. For the constant £ case (right) the multiphase strength is given by eg (equation ^) and 
for (left) the 5-function case cr? (equation M) for N = 2 and N c = 20. 



erg which is just AT/2{T) for a constant £. Again, for the Si 
ratio there is a (T) near 1 keV where a large range of <rg have 
nearly the same Hl/He4 ratio; i.e. even for this rigid model, 
more than one line ratio is required to constrain eg . Notice, 
however, that there are large regions of ratio- (T) space which 
are uninhabited. In particular, the accessible region of the 
£ = constant model ends not much below the line erg = 0.9. 
This also applies for any continuous £(T) since £ = constant 
is defined to be the limiting case for maximal multiphase 
strength. (Note that a Gaussian model has behavior similar 
to that of the constant £ in Figure ^) 

From examination of plots like Figure ^ for different 
elements we can locate the critical value of ag that rep- 
resents the boundary between approximately single-phase 
and strongly multiphase £. If we define £ to be significantly 
multiphase when the Hl/He4 line ratios of Si (and others) 
typically differ by more than 10 per cent from the isother- 
mal case, then erg ~ 0.2 represents the boundary between 
approximately single-phase and significantly multiphase £. 
This criteria applies for a constant and Gaussian £. When 
a 5-function is added to these models we find that erg ~ 0.1 
is a better representation of the transition between a single- 
phase and multiphase gas. 

2.2.3 Cooling flows 

A physical mechanism for producing multiphase gas in el- 
lipticals, groups, and clusters is provided by cooling flows 
(Fabian, Nulsen, & Canizares 1984; Nulsen 1986, 1997; 
Thomas, Fabian, & Nulsen 1987; White & Sarazin 1987; 
for a recent review see Fabian 1994). In these models gas 
cools out of the flow in regions where the cooling time 
is sufficiently short (usually taken to be less than the as- 
sumed age of the system). The simplest example is for gas 
cooling at constant pressure from some upper temperature 
Tmax in which case the conservation of enthalpy implies 
an energy release dL within a temperature interval dT of, 
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Figure 7. The differential emission measures for a constant pres- 
sure cooling flow (Z = O.5Z0 - solid, Z = 2Zq — long dash) and 
for a cooling flow with M oc r in the gravitational potential of 
a singular isothermal sphere (Z = O.^Zq — dotted, Z = 2Zq - 
short dash). All models have T max = 4 keV. The units of £ (for 
each metallicity) are normalized to the maximum value of £ for 
the constant pressure cooling flow. 

dL = 5A/fcfldT/2/im p , where /i is the mean atomic weight, 
m p the proton mass, and M is the rate of mass deposition. 
Comparing this expression to equation (^) one obtains the 
differential emission measure, 

5Mfcs_l (Constant Pressure CF), (8) 
2 /!% A(T) 

for T < T m ax which we plot in Figure for T max = 4 keV 
with abundances Z — 0.5, 2Zq. The shape of f (T) is sensi- 
tive to Z, particularly for T ^ 1 keV where the Fe L shell 
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Figure 8. Multiphase strengths for the cooling flow models in 
Figure [?j. See text for discussion of the behavior near 1 keV. 



lines dominate A(T). However, since most of the emission 
measure, J £(T)dT, in the models in Figure ^ result from 
temperatures between ~ 2 - 4 keV, these metallicity differ- 
ences do not translate to large differences in the K shell line 
ratios as we show below. 

The emission measure weighted temperatures for the 
models with Z = 0.5/7q range from (T) ~ 0.6T max for 
T max = 0.5 keV to (T) « 0.5T max for T max = 10 keV. The 
luminosity weighted temperatures are w 15 per cent larger 
than (T). For Z = 2Z Q , (T) ranges from (T) ^ 0.4T max for 
= 0.5 keV to (T) w 0.8T max for T max = 10 keV. 

The multiphase strength as a function of T max is shown 
in Figure |§ Over T max = 0.5 - 10 keV we have a € ~ 0.45 
- 0.8 for both Z — 0.5Zq and Z = 2Zq and thus the con- 
stant pressure cooling flow has nearly maximal multiphase 
strength for the temperatures of ellipticals, groups, and clus- 
ters. The Z — 2Zq case lies below the Z = 0.5-Zq model but 
they approach each other at the largest values of T max . Note 
that the wiggles in <j£ are simply due to the same wiggles in 
A. 

In Figure we show the Hl/He4 line ratios of the con- 
stant pressure cooling flow (with Z = 0.5Zq) as a function 
of T max . The Hl/He4 profiles are markedly different from 
the single temperature case shown in Figure |§| The differ- 
ences are most striking for those lines with E ^ 2 keV. In 
particular, the O Hl/He4 profile changes by only a factor 
of ~ 3 while the single temperature profile changes by sev- 
eral orders of magnitude over the temperature range shown. 
(The O lines are much stronger in the cooling flow models 
as well.) In addition to the K shell lines, the Fe L line ra- 
tios (not shown) also have very different profiles from the 
isothermal case making them very powerful for constraining 
the properties of cooling flows (e.g. Canizares et al. 1982, 
1988). 

The line ratios are very similar for the Z = 2Zq case 
where the Ca and Fe Hl/He4 profiles deviate from those 
in Figure [] by less than ~ 5 per cent over the whole tem- 
perature range plotted. The O and Si ratios show larger 
deviations which increase with increasing T max ; i.e. above 2 



keV the O and Si ratios are larger by ~ 20 - 40 per cent and 
~ 10 - 20 per cent respectively for the Z = 2Zq case. 

Since only regions where the cooling time is less than 
the age of the system participate in the cooling flow, the 
observed X-ray flux will generally include emission from the 
non-cooling "ambient" phase. In Figure ^| over the constant 
pressure model we also plot the Hl/He4 ratios of a model 
consisting of a constant pressure cooling flow of tempera- 
ture, Tmax and a single temperature model with tempera- 
ture, T = T max where the cooling-flow and isothermal com- 
ponents have equal total emission measure. The line ratio 
profiles of the composite models have reasonably similar- 
shapes to the case of the cooling flow alone, but they are 
shifted up by ~ 10 - 40 per cent. The largest differences oc- 
cur at the largest temperatures for the lines with the highest 
energies. 

The reason why the composite model still has line ra- 
tios markedly different from the isothermal case is that each 
component has the same total emission measure which im- 
plies / = 0.5 and Of — > 0.5<7£ for the co mposite cooling flow 
and S- function model (see section [2.2.1 ); i.e. the values of crj 
range from ~ 0.3 - 0.4. As noted at the end of the previous 
section, the boundary between approximately single-phase 
and significantly multiphase emission measure distributions 
occurs for <j£ ~ 0.1 when £ includes a J- function compo- 
nent. Thus, a model with a cooling flow component and a 
single-phase component has line ratios similar to an isother- 
mal gas for / <, 1 /7 and E1M5 /EM cf <; 6. (Note this applies 
for Z = 0.5Zq but is similar for Z — 2Zq where we deduce 
/ <, 1/5 and EM 4 /EM cf J> 4.) 

It is worth considering how much these properties 
change for a more sophisticated multiphase cooling flow 
model that allows for g ravita tional work to be done on the 
flow as it cools. Nulsen (1998) has derived £(T) for a cooling 
flow in the gravitational potential of a singular isothermal 
sphere. For T <C Tm ax the differential emission measure of 
this isothermal sphere cooling flow model is equivalent to the 
constant pressure case (equation]^). In contrast, £(T max ) = 
for the isothermal sphere cooling flow. The detailed shape 
for intermediate temperatures depends on the slope rj of the 
mass deposition profile, M « r'. 

We plot £(T) and for rj = 1 in Figures |^ and ^| for 
comparison to the constant pressure case. Except for T max ~ 
1-2 keV, the multiphase strength for the Z = 0.5Zq case 
is slightly smaller for the isothermal sphere cooling flow; 
the differences are more drastic for the Z = 2Zq case. The 
emission measure weighted temperatures for these models 
are ~ 5 - 10 per cent smaller than those of the constant 
pressure cooling flow. 

There is an abrupt change in er^ near 1 keV for the 
isothermal sphere cooling flow models. Let us focus on the 
Z — 0.5Zq case first. Near 1.2 keV £ max jumps from its 
location near the local maximum at 0.36 keV to a tempera- 
ture just below T max (see Figure . As T max increases over 
the next ~ 0.3 keV we have £ max increasing faster than the 
total emission measure translating to a decrease in <7{. This 
change is much more abrupt for the isothermal sphere cool- 
ing flow because of the enforced condition that £(T max ) = 0. 
This overall behavior is more pronounced for the Z = 2Zq 
case because of the details of the strong Fe L shell lines are 
proportionally more important. 
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Figure 9. Hl/He4 ratios for selected elements in models with a constant pressure cooling flow (left) and isothermal sphere cooling 
flow with M at r (right). In each case the line ratios for the cooling flows are given by solid lines while composite models consisting of 
a cooling flow and a single temperature component (each with equal total emission measure) are given by the dotted lines. All models 
have Z = 0.5Zq and relative abundances are fixed at the (photospheric) solar values. 



The Hl/He4 line ratios for both the isothermal sphere 
cooling flow component and a composite model of an isother- 
mal sphere cooling flow component with a single tempera- 
ture component (all with Z = 0.5Zq) are plotted along 
side the corresponding constant pressure models in Figure 
^| For the case of only a cooling flow the line ratios are very 
similar, especially for T ma x is 2 keV. The greatest discrep- 
ancies between the isothermal and constant-pressure cases 
are seen at the largest temperatures for the lines with the 
highest energies; e.g. the Fe Hl/He4 ratios differ by ~ 40 
per cent for T max = 10 keV. Adding the single temperature 
to the isothermal cooling flow has a slightly larger effect 
than found for the constant pressure cooling flow. However, 
the composite models of the isothermal and constant pres- 
sure cases actually differ from each other by ^ 10 per cent 
over the entire temperature range shown. (Note that the 
line ratios for models with Z = 2Zq differ by essentially 
the same small amounts described above for the constant 
pressure models.) 



2.3 Distinguishing a two-temperature plasma 
from a cooling flow 

We have seen in the previous section that for a particular 
value of a line ratio (e.g. Si Hl/He4) the two-temperature 
models appear to have a larger variety of possible which 
can produce a line ratio in comparison to the simple con- 
tinuous emission measure models examined. As a result, for 



moderate resolution spectra like ASCA ( Fanaka et al. 1994 ) 



where at most only a limited number of line blends are re- 
solved, it should be expected that a two-temperature model 
can provide an acceptable fit even if the true spectrum is 
that of a cooling flow. In fact, analyses of ASCA spectra of 
the brightest e llipticals (Buo te & Fabian 1998; Buote 1999a), 
poor groups (Buote 19991), and clusters of galaxies (e.g. 
Fabian et al. 1994a) have found that two-temperature mod- 
els and models composed of a cooling flow with an ambient 



temperature component essentially describe the data equally 
well. 

To examine how to distinguish between the two- 
temperature and cooling flow models we first examine the 
case of the bright elliptical galaxy NGC 4472. In Table ^ we 
list the results of the best-fitting two-temperature and cool- 
ing flow models to theASCA spectra accumulated within 
a radiu s of ~ 30 kpc from the center of NGC 4472 by 
Buote (p.999a|). The "2T+BREM" model consists of two 



MEKAL plasma components and one component of ther- 
mal bremsstrahlung to account for possible emission from 
discrete sources; T c and Th represent the temperatures of 
the "cooler" and "hotter" MEKAL components respectively. 
The "CF+1T+BREM" model replaces the cooler MEKAL 
component with a constant pressure cooling flow component 
(equation |^). Each model component is modified by photo - 
electric absorption (Balucihska- Church & McCammon 1992) 



although the absorption on the bremsstrahlung components 
is fixed to the Galactic value. We refer the reader to Buote 



(1999a) for full details of these models. 

In Figure [n] we plot a selection of line ratios computed 
from the models in Table (We show a related plot of the 
ratios of the spectra of two temperature and cooling flow 
models in Figure ^ in Section [l^.) Overall, it is remark- 
able how similar the line ratios are for the two models. The 
O Hl/He4 and Fe L ratios show the most differences be- 
tween the two-temperature and cooling flow models; i.e. the 
O Hl/He4 ratios differ by a factor of ~ 4 and the 2:1 Fe 
L ratios differ by a factor of ~ 2. The ratios shown involv- 
ing other elements differ by ^ 10 per cent. Note that we 
expect the Hl/He4 ratio of Si to be similar for the models 
since the blends involving the He4 and HI transitions are 
reasonably resolved by ASCA (see next section) and thus 
the best-fitting models in Table attempt to match the ob- 
served Si ratio as accurately as possible. 

For comparison to NGC 4472 where the cooling flow 
component dominates the emission measure and T max < 2 
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Table 4. Spectral Models for NGC 4472 and the Centaurus Cluster 

7V£ T c T h Z EM C EM h M T B EM B <r c 

Model (10 21 cm- 2 ) (keV) (keV) (Z Q ) (10 65 cm- 3 ) (M Q yr" 1 ) (keV) (10 65 cm" 3 ) 

NGC 4472: 

2T+BREM 2.9 0.04 0.72 1.39 2.00 0.065 0.050 ■•■ 3.23 0.0069 0.44 

CF+1T+BREM 2.9 0.00 1.30 1.30 1.41 0.185 0.032 2.54 3.71 0.0045 0.03 

Centaurus: 

2T 2.3 1.4 1.04 3.07 1.10 1.41 17.0 ■•■ ■■• ■■• 0.02 

CF+1T 3.0 1.3 3.30 3.30 1.19 2.45 11.8 45.4 ■■• ■■• 0.11 

Best-fitting spectral models to ASCA data for the elliptical galaxy NGC 4472 and for the Centaurus cluster. The CF components are 
constant pressure cooling flow models, the IT components are single-temperature MEKAL models, and B stands for a bremsstrahlung 
component. See text for further explanation of these models. The emission measures are computed using the distances described in the 
notes to Table |. (Note that ct 5 for the 2T and 2T+BREM models is actually evaluated for N c = 20.) 
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Figure 10. Selected line ratios of two-temperature models (boxes) and models consisting of a cooling flow and a single temperature 
(crosses) for NGC 4472 and Centaurus (see Table [ij). (Note that the line ratios are computed as described in Section ^ except that the 
absorption is included.) 



keV we consider the Centaurus cluster of galaxies. In Ta- 
ble we list the best-fitting 2T and CF+1T models to the 
ASCA data accumulated within a radius of ~ 100 kpc. The 
details of the data reduction and spectral fitting follow those 
in section | with the following important exceptions. Firstly, 
we fitted the ASCA SIS spectra over the energy range 0.5-9. 
keV. Secondly, we used the phot osphe ric solar abundances 
according to Anders & Grevesse (1989) because of our cur- 
rent im plementation of the cooling flow model for use in 
XSPEC (Arnaud 1996) is defined for those solar abundances. 



Although this restriction limits the quality of the fits, the 
2T and CF+IT models give very similar quality fits which 
is all that we require for this discussion. We refer the reader 
to section |^ for a thorough discussion of successful spectral 
models of Centaurus. 

In contrast to NGC 4472 the isothermal component 
dominates the emission measure of the cooling flow model 
for Centaurus, EM h /EM c » 5. Nevertheless, the CF+IT 
model has crj = 0.11 indicating that it is significantly multi- 
phase (see previous section). We plot selected line ratios of 
the 2T and CF+IT models in Figure [lo| Since Centaurus 
has a higher temperature than NGC 4472, in the figure we 



replace the S Hl/He4 and Sat/He4 ratios with the corre- 
sponding Fe ratios; S/Si (He4) is replaced with S/Fe (He4). 
(Again, we show a related plot in Figure [l^ in Section 4.2.) 
Similar to NGC 4472 the O Hl/He4 and Fe L ratios show 
the most substantial differences (factors of ~ 5). However, 
unlike NGC 4472 the O He3/He4 and Si Sat/He4 show no- 
ticeable differences. (The Sat/He4 ratios of S, Ar, and Ca 
have quite similar differences as well.) Again the Si Hl/He4 
ratios are very similar for the three models because the Si 
Hl/He4 ratio is a key constraint of the ASCA spectral fits. 

Hence, the most substantial differences between two- 
temperature and cooling flow models occur for O and Fe 
L line ratios. Analyses that exclude energies below 1.4 keV 
(e.g. to avoid small inaccuracies in the plasma codes) will 
have a very difficult time differentiating two-temperature 
and cooling flow models, particularly for moderate spectral 
resolution data provided by ASCA, Chandra, and XMM. 
Moreover, because of their overall flexibility, an observed 
consistency of a two-temperature model with data of low- 
moderate spectral resolution of ellipticals, groups, and clus- 
ters should not be regarded as decisive evidence that the 
system has only two phases. 
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Table 5. Energy Ranges of ASCA Line Blends 

Energy Range 
Dominant Ion (keV) 

Si XIII 1.80 - 1.90 

Si XIV 1.95 - 2.05 

S XV 2.35 - 2.53 

S XVI 2.55 - 2.70 

Ar XVII 3.00 - 3.25 

Ar XVIII 3.25 - 3.40 

Fe XXV 6.40 - 6.80 

Fe XXVI 6.86 - 7.14 

7.8 keV Blend 7.60-8.10 

Energy ranges used to define the lines in Figure [ll] appropriate 
for the energy resolution of the ASCA SIS during approximately 
1994. 



3 APPLICATION TO ASCA DATA 

Recently the ASCA satellite is able to resolve emission line 
blends of K shell transitions of Mg, Si, S, Ar, Ca, and Fe, and 
also has the effective area and imaging capabilities lacking 
in previous X-ray satellites having similar or better spectral 
resolution. The energy resolution of ASCA (e.g. AE ~ 100 
eV FWHM at 2 keV) allows only blends of lines to be mea- 
sured, and we show in Figure [ll] the temperature profiles of 
ratios of the most important line blends relevant to ASCA 
observations of ellipticals, groups, and clusters; the defini- 
tions of the blends are given in Table []. These profiles should 
be compared to Figure ^[ 

Even though they are the most p owe rful diagnostics of 
multiphase cooling flows (see Section 2.3), we do not show 
oxygen or Fe L blends in Figure [ll] for the following reasons. 
The strongest transitions of the He-like O VII ion appear at 
energies ~ 0.57 keV which lie in a region where the ASCA 
SIS have small effective area and problems with the calibra- 
tion. The energy resolution of ASCA does not allow useful 
separation of blends of Fe L lines between ~ 0.7 - 1.4 keV for 
individual analysis. This is not crucial to our present analy- 
sis since it is our intention to investigate the K shell emission 
lines which do not suffer from the possible inaccuracies in 
the plasma codes associated with the Fe L shell transitions. 

As indicated by the discussion in Section 2.Z line ratios 
other than those of oxygen or Fe L in rival multiphase models 
generally differ by only ~ 10 - 20 per cent. Only a handful of 
ASCA observations of the brightest ellipticals, groups, and 
clusters have sufficient S/N to even approach this level of 
precision. The Ka line blends of both H-like and He- lik e Si 
and S have been analyzed with ASCA by Buote (1999a) for 



three of the brightest elliptical galaxies (NGC 1399, 4472, 
and 4636) and one of the brightest galaxy groups (NGC 
5044). Although the properties of these line blends favor 
multitemperature models obtained from broad-band (0.5- 
9 keV) fits to the ASCA spectra, the line blends are not 
measured precisely enough to place interesting constraints 
on £(T) by themselves. 

Clusters of galaxies are hotter and more luminous 
than ellipticals and groups, and thus the K shell emis- 
sion lines above ~ 2 keV in principle can be measured 
more precisely for the brightest clusters. Let us consider 
the bright nearby clusters M87 (Virgo), Centaurus (A3526), 
and Perseus (A426) which have among the highest S/N 
data available for clusters observed with ASCA. Broad-band 



spectral fitting of the ASCA data (which includes the Fe L 
lines) has shown that at least two temperature components 
are required in the cores of each of these clusters (Fabian et 
al. 1994a; Fukazawaet al. 1994; Matsumoto et al. 1996). It is 
our intention to examine what can be learned from analysis 
of only the K shell emission line blends and, in particular, 
to examine whether they too require multiple temperature 
components in M87, Centaurus, and Perseus. 



3.1 ASCA observations and data reduction 

The ASCA satellite consists of two X-ray CCD cameras 
(Solid State Imaging Spectrometers - SIS0 and SIS1) and 
two proportional counters (Gas Imaging Spectrometers - 
GIS2 and GIS3) each of which is illuminated by its own 
X-ray telescope (XRT). The SIS has superior energy reso- 
lution and effective area below ~ 7 keV while the GIS has 
a larger field of view. Although the point spread function 
(PSF) of each XRT has a relatively sharp core, the wings of 
the PSF are quite broad (half power diameter ~ 3') and in- 
crease markedly for energies above a few keV (e.g. Kunieda 
et al. 1995). As in our previous studies we do not attempt to 
analyze the spatial distribution of the ASCA data of these 
sources because of the large, asymmetric, energy-dependent 
PSF. Rather, we analyze the ASCA X-ray emission within 
a single large aperture for each cluster which encloses the 
region of the brightest emission (see below). For clusters of 
galaxies the spectrum within this central region is not very 
much affected by scattered emission from large radii (e.g. 
Boehringer et al. 1998). 

We obtained ASCA data for the clusters from the public 
data archive maintained by the High Energy Astrophysics 
Science Archive Research Center (HEASARC). The prop- 
erties of the observations are listed in Tables ll and pj. We 
reduced these data with the standard FTOOLS (v4.1) soft- 
ware according to the procedures described in The ASCA 
Data Reduction Guide and the WWW pages of the ASCA 
Guest Observer Facility (GOFQ 

For the observations of M87 and Perseus (taken during 
the Performance- Verification phase) we used the events files 
generated by the default screening criteria processed under 
the Revision 2 Data Processing (REV2) . As is standard pro- 
cedure for analysis of SIS data, only data in BRIGHT mode 
taken in medium or high bit rate were used. Since the Cen- 
taurus observation was performed in 1995 its SIS data are 
degraded to some extent because of radiation damage to 
the CCDs (see Dotani et al. 1997). We corrected the Cen- 
taurus data for this "RDD effect" following the procedure 
outlined on the ASCA GOF WWW pages which produces 
corrected SIS event files in BRIGHT2 mode constructed us- 
ing the default REV2 screening criteria. (Again, only data 
with medium or high bit rate were used.) Since we focus on 
energies above 1.6 keV this correction is actually not very 
important for our study. 

Our final screening of the data involved visual exami- 
nation of the light curves for each observation and removing 
intervals of anomalously high (or low) count rate. We note 
that only data with the standard events grades (0234) were 



See |http: / /hcasarc. tJsfc.nasa.gov/docs/asca/abc/abc. html 



and rttp:/ /heasarc. gsfc.nasa.gov/docs/asca/ 
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Figure 11. (Left) Line ratios as a function of temperature in isothermal MEKAL models (with Z = 0.5Zq) where the emission lines 
have been computed over large energy ranges appropriate for the energy resolution of ASCA (see Table The ratios for Si, S, Ar, and 
Fe refer to ratios of H-like to He-like Ka transitions as in Figure ^ except that now additional transitions are included in the blends. The 
7.8 keV blend consists mostly of K/3 Fe XXV (He3) although Ka emission from Ni XXVII and Ni XXVIII contribute as well. (Right) 
The same thing but plotted on a linear scale. 



Table 6. Cluster Properties and ASCA Observations 



Cluster 


z 


N n 


Sequence 


Date 


SIS 


SIS 


RDD 


Radius 


Radius 






(10 21 cm~ 2 ) 




Mo/yr 


CCD 


Data 


Correction 


(arcmin) 


(kpc) 












Mode 


Mode 




SIS 


GIS 


SIS GIS 


M87 


0.0043 


0.25 


60033000 


6/93 


4 


BRIGHT 


N 


6.2 


7.5 


32.5 39.3 


Centaurus 


0.0110 


0.81 


83026000 


7/95 


1 


BRIGHT2 


Y 


5.6 


6.5 


77.2 89.6 


Perseus 


0.0183 


1.41 


80007000 


8/93 


1 


BRIGHT 


N 


5.6 


6.3 


128.7 144.8 



Select ed properties of the ASCA SIS and GIS observations. Galactic Hydrogen column densities (Nn) are from Dickey & Lockman 
(1990) using the HEASARC w3nh tool. Extraction radii in kpc are computed using luminosity distances of 18 Mpc for M87, 47 Mpc 
for Centaurus, and 79 Mpc for Perseus assuming H = 70 km s" 1 Mpc -1 and fl = 0.3. 



used, and we made the required dead time corrections for 
the GIS data. 

The final processed events were then extracted from a 
region centered on the emission peak for each detector of 
each sequence. We selected a particular extraction region 
using the following general guidelines. Our primary concern 
is to select a region that encloses most of the X-ray emission 
yet is symmetrically distributed about the origin of the re- 
gion. We limited the size of the aperture to ensure that the 
entire aperture fit on the detector, an issue more important 
for the SIS because of its smaller field-of-view. Moreover, 
for the SIS0 and SIS1 we tried to limit the apertures to as 
small a number of chips as possible to reduce the effects 
of residual errors in chip-to-chip calibration. The GIS re- 
gions were chosen to be of similar size to the corresponding 
SIS regions of a given sequence for consistency. Since, how- 
ever, the GIS+XRT PSF is somewhat larger than that of 
the SIS+XRT the extraction regions for the GIS are usually 
~ 20% larger. 

Since the SIS was operated in 1-ccd mode for the Cen- 
taurus observation we extracted all of the emission from that 
ccd which amounts to a square region with sides of ~ 11' 
width. The observation of Perseus was positioned very near 
the center of the default chips of the SIS and so we also 
extracted its SIS data from the whole chip. The observa- 



tion of M87 was centered between two chips, and thus we 
used a circular region. In Table ^| we list the region sizes for 
the SIS and GIS data of each cluster. For the SIS data of 
Centaurus and Perseus the radii actually are half-widths of 
the default ccds. (Note that the data between the ccds are 
excluded.) All of the GIS radii refer to circles. We extracted 
the events using regions defined in detector coordinates as is 
recommended in the ASCA ABC GUIDE because the spec- 
tral response depends on the location on the detector not 
the position on the sky. 

We computed a background spectrum for each detec- 
tor using the standard deep observations of blank fields. It 
is necessary to use these background templates because the 
cluster emission fills the entire field of view in each case The 
background templates for the SIS are most appropriate for 
SIS data taken in 4-ccd mode early in the mission which 
applies to the M87 and Perseus observations. For more re- 
cent data taken in 1-cdd or 2-ccd mode, using these back- 
ground templates can lead to spurious effects particularly 
for energies above ~ 7 keV where instrumental background 
dominates the cosmic X-ray background. The Centaurus ob- 
servation was performed in 1995 and thus these effects are 
not expected to be as serious as for more recent data. At any 
rate, the Centaurus data are not dominated by background 
so small systematic effects in the background should not be 
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Table 7. ASCA Exposures and Count Rates 



Cluster 


Exposure 


Count Rate 


Exposure 


Count Rate 




(10 3 


s) 


(ct s" 1 ) 


(10 3 s) 


(ct s- 1 ) 




SISO 


SISl 


SISO SISl 


GIS2 GIS3 


GIS2 GIS3 


M87 


12.8 


9.79 


2.59 1.87 


16.7 16.7 


2.33 2.39 


Centaurus 


64.6 


64.4 


1.22 0.94 


70.8 70.8 


0.86 1.06 


Perseus 


16.1 


16.1 


7.15 6.16 


15.0 14.9 


5.53 6.62 



Exposures include any time filtering. The 1.6-9 keV count rates are background-subtracted within the particular aperture (see Table o 



apparent in our analysis. The final screened background- 
subtracted exposures and count rates are listed in Table [?]. 

The instrument response matrix required for spectral 
analysis of ASCA data is the product of a spectral Redistri- 
bution Matrix File (RMF) and an Auxiliary Response File 
(ARF). The RMF specifies the channel probability distribu- 
tion for a photon (i.e. energy resolution information) while 
the ARF contains the information on the effective area. An 
RMF needs to be generated specifically for each SIS of each 
observation because, among other reasons, each chip of each 
SIS requires its own RMF and the spectral resolution of the 
SIS is degrading with time. We generated the responses for 
each SIS using the FTOOL SISRMG selecting for the stan- 
dard event grades (0234). Using the response matrix and 
spectral PI (Pulse Invariant) file we constructed an ARF 
file with the FTOOL ascaarf. 

The source aperture used for the SIS overlapped more 
than one chip for M87. In order to analyze the spectra of 
such regions we followed the standard procedure of creating 
a new response matrix that is the average of the individual 
response matrices of each chip weighted by the number of 
source counts of each chip (i.e. within the source aperture). 
(Actually, the current software only allows the RMFs to be 
averaged and then an ARF is generated using the averaged 
RMF. This is considered to be a good approximation for 
most sources.) Unfortunately, some energy resolution is lost 
as a result of this averaging. 

Since the K shell emission lines for energies above 
~ 2 keV are relatively faint in these clusters we maximize 
the S/N of the data by combining the spectra from each 
SIS and similarly for each CIS; i.e. SIS=SIS0+SIS1 and 
GIS=GIS2+GIS3. For the CIS data this does not result 
in significant loss of information because the RMFs of the 
GIS2 and GIS3 are the same and time-independent; only the 
ARFs need to be averaged for the GIS2 and GIS3. However, 
the responses of the SISO and SISl detectors are different 
and thus the required averaging of the RMFs and ARFs 
leads again to loss of information. Because the gain in S/N 
is important for our study, we examine the combined SIS and 
combined GIS data of each cluster where the source spectra, 
RMF, ARF, and background for each detector are summed 
and scaled appropriately using the FTOOL ADDASCASPEC 
(vl.27). 

Overall we found that the results obtained from analysis 
of these summed SIS spectra agreed very well with results 
obtained without summing the data. Henceforth we shall fo- 
cus on results obtained from the summed SIS data though 
we will (for completeness) mention some of the results ob- 
tained from separate analysis of the SISO and SISl data. 



3.2 Analysis of individual emission line blends 

We focus our analysis of individual line blends on the SIS 
data because, unlike the GIS, the SIS is able to resolve both 
the He-like and H-like blends of the K shell lines of the 
abundant elements (e.g. see Figure 1 of Buote 1999a). The 
GIS data can provide useful supplementary constraints on 
the continuum for energies above ~ 5 keV which we mention 
below. 

Our procedure for obtaining the properties of the line 
blends commences by fitting a thermal bremsstrahlung 
model to represent the continuum emission. We desire a 
measurement of the local continuum near the line blends 
of interest so that the assumption of a single temperature is 
reasonable even if the spectrum actually consists of multi- 
ple temperature components. However, because of S/N con- 
siderations and the finite energy resolution of the SIS we 
allowed the continuum component to extend over multiple 
line blends that are nearby in energy. The line blends of 
interest were then modeled as Gaussians of zero intrinsic 
width on top of the continuum component. 

We fitted only line blends that are either obvious from 
visual inspection of a spectrum or significantly improved 
the fit as judged in terms of \ 2 . All spectral fitting was per - 
formed with the software package XSPEC (Arnaud 1996). 
In order for the weights to be valid for the x z method we 
followed the standard procedure and regrouped the PI bins 
for each source spectrum so that each group had at least 
20 counts. The background templates also generally have at 
least 20 counts when their energy bins are grouped similarly 
to the source spectra. The background-subtracted count rate 
in a particular group can be small, but the uncertainties 
in the source and background are correctly propagated by 
XSPEC to guarantee approximately gaussian statistics for 
the statistical weights (K. Arnaud 1998, private communi- 
cation). 

The energies of the lines were free parameters in the fits. 
When determining the error bars on the equivalent widths 
and fluxes for a particular line the energies of the other 
lines were fixed at their best-fitting values. This was done 
to insure that the other lines did not mix during the error 
search. 



3.2.1 Properties of Ka and Kf3 lines 

In Table ^| we summarize the properties of the line blends 
measured from the summed SIS data. The spectral regions 
and best-fitting models used to measure the line properties 
are displayed for each cluster in Figures |il||l4| . 

We detect most of the strong Ka lines of the abun- 
dant elements lying between the energies 1.6 and 9 keV in 
each cluster . In Centaurus both the H- and He-like Ka line 
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Table 8. Emission Lines Measured from ASCA Data 

E T EW Flux 

Ion (keV) (keV) (eV) (1CT 5 ph cm" 2 s" 1 ) (1CT 13 erg cm" 2 s" 1 ) 

M87: 

MgXII(H) 1.48±S:8a l-06±g;i| 28±| 62.86±£°£ 14.86j£g 

1.7 keV Feature 1.74+°-™ 1.23+°-^ 13±| 19.91±|"^ 5 - 54 i255 
Si XIII (He) 1.84±g : gJ 1.23±g : » 54±g 68.68^ : g 20.29±^° 
Si XIV (H) 1.99lg;gg 1.23lg:g 110±1J U0.4D±&gg 36.26l|;g 
2.2 keV Feature 2.18±g-jg 1.23±g-£ 28±E 21.94l 10 97 7 7.66±|-g 
SXV(He) 2.44±g : g} 2.19l° : ^ 68±i* 47-25j*£ 18.49tjj;°J 
S XVI (H) 2.60!°'°^ 2 - 19 ±o'bo 68 -n 40.52±HI 16.90±|'^ 
Ar XVII (He) 3.Q9±g;|g 2.43l°;^ 2 44^ 2 16.86±|;|| S- 35 ^!? 
*Ar XVIII (H) 3.28l°,;°J 2.43±g;^ 111 7 , 3.541?,;^ i-^tiiw 
Ca XIX (He) 3.89lg : g| 2AZ+°f 9 43l^ 8.70±|*|| 5.42±^ : |g 
FeXXV(He) 6.63±g;g* 2.61±g;g 1010±E° 32.57±t?? 34.61±jp 

7.8 keV Feature 7.82±g;^j| 2.61±g;g 308jjg! 5.02j;|;|| 6.29l*;°° 

Centaurus: 

MgXII(H) 1.47±g-jg 3.06lJ-°2 2711 26.22±|'|5 6.18l°-™ 

Si XIII (He) 1.81±g gi 1.6llg ^ 2511 16.521^° 4.781°,^ 

Si XIV (H) 1.98l°- °° 1.6ll°'jg 83l* 42.431 3 13.47l{'^ 

+0:02 +o'.24 S4+ 13 13 17+ 4 ' 26 4 68 +1 ' Sl 



2.2 keV Feature 2.22l°'g 2 1.61±g-J| 34±^| 13.171* 4.681^^ 

SXV(He) 2.44l° : °i 2.97l o : °? 44l g ° 16J}2±| : gl 6.4B±li» 

S XVI (H) 2.601°/°^ 2.97lJ'°* 501* 16.09l|'^ 6 - 71 io 90 

Ar XVII (He) 3.09l^ 4.12+g* 40i° 8.58±};£ 4.25+_° f 9 

Ar XVIII (H) 3.31±g-{g 4.12±g-£ 22l* 4.15±°^| 2.20l°'™ 

CaXIX(He) 3.86±° : ^ 4.12l° : « 26l^ 3.54±*J° 2.191 ,-' 7 , 1 

CaXX(H) 4.04lg;g? A.12±H 12l 9 8 1.53±J|g| 0.99!°;™ 

Fe XXV (He) 6.63l°;°i 4.37±g;|| 972lg 33.8ll*j 3 35.891^ 

Fe XXVI (H) 6.971°;°^ 4.371°;^ 4ll 2 * 2 .40l^ 2.68+^;|f 

7.8 keV Feature 7.8ljlg;^ 4.37±g;|| 211±|* 4.21±^ 5.26±^|? 

8.4 keV Feature 8.35lg;°® 4.37lg;^ 24l!J 2 ® 3.86l^™ S.lTl^ls 

74 1Q+16.28 i 7fi 2+ 3 ' 87 

'*- la -lS.07 A '* DZ -3.58 

42 i , i+ 18 ' 22 ii fin+ 5 02 

^• 1J _18.ao ii ' OU -5.01 

54 45+61-35 ioq2+ 3 ' 38 

101.73t»:« 32.18l£S 

75.19±J|-|| 26.82lS-^ 

47.15l 12 i ;°« 19.68±|i°* 

1S2^ +9,53 6 00+ 432 

10.ZO_g 30 D.UU_ 4 21 

21.08±|-.|| 10.451^ 

IS.61!^ 1 9.72l^ 3 

UA0+_tf 7 6.971 3 ™ 

9.32l^ 6.03l«;« 

157.71±|-g 166-S5+|-J| 

37.14±™ 4D.54±|;iS 

13.55t*;gJ 16.80lg : » 

Results of fitting thermal bremsstrahlung and Gaussian models to ASCA SIS data of ellipticals and clusters having the largest available 
S/N. Except for Fe XXV (He) the intrinsic widths of the Gaussians are set to zero. Lines that have the same continuum temperature, 
T, are cases where multiple gaussians are joined by one bremsstrahlung model. See text for additional details regarding the model 
fitting. Quoted errors are 90 per cent confidence on one interesting parameter except for those lines with an asterisk for which 68 per 
cent limits are quoted. 
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blends of Si, S, Ar, Ca, and Fe are detected. These lines are 
also detected in the other clusters except for S XV (He) in 
Perseus and for Ca XX (H) and Fe XXVI (H) in M87. Mg 
XII (H) is also detected in each case but we advise caution 
in interpreting the derived EWs and fluxes of this line be- 
cause the nearby unresolved Fe L lines undoubtedly affect 



our estimates of the local continuum and the flux arising 
from only Mg. 

In addition to the expected strong Ka lines we detect 
other significant lines in each cluster. We label these as "Fea- 
tures" in Table [] because they tend to be weaker than the 
Ka lines and, more importantly, some of these line features 
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Figure 12. ASCA SIS spectral data for M87. The top four panels sh ow the energy ranges and best-fitting Gaussian+continuum models 
used to obtain the line properties listed in Table ^ See text in section ^2] for further explanation of these models. The bottom two panels 
show the best-fitting isothermal (left) and two- te mpe rature (right) models fitted to the SIS data over 1.6-9 keV as listed in Table [lo[ 
These broad-band models are discussed in section 3.3. 



may result in part or entirely from systematic errors in the 
SIS response. In particular, for each cluster we measure a sig- 
nificant line feature near 2.2 keV which in principle could be 
attributed to the 2.18 keV K/3 line Si XIII (He3). However, 
this is precisely the location of a known calib ration prob- 
lem ow ing to the optical constants in the XRT (Gendreau & 
Yaqoot 1997) making it difficult to determine reliably the 



origin of the 2.2 keV features. 



A significant line feature is also detected near 1.7 keV 
in M87 and Perseus. Assuming the feature is not instrumen- 
tal then it is a blend of KB Mg (H), Ka Al (H), and some 
Fe L emission (see Figure pi). Although we are not aware of 
any problems currently in the SIS response near 1.7 keV, 
we stress that these features (and those at 2.2 keV) are not 
especially prominent (see Figures [l^ and |l4j) and their exis- 
tence is thus very dependent on the calibration. 
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Figure 13. As Figure |l^ but for Centaurus 



In each cluster we also detect line features near 7.8 keV. 
If these features originate from hot plasma, they are primar- 
ily the result of K/3 emission from Fe XXV (He3) although 
some Ni Ka should contribute as well. This feature is most 
prominent in Centaurus though it is clearly apparent to the 
eye in the SIS spectra of M87 (though it is noisy) and is 
least noticeable for Perseus (see below). Visual examination 
of the M87 spectrum (Figure |l^) near 8 keV suggests that 
the line feature is broader than the instrumental resolution. 
A slightly better fit can be achieved either be adding an- 
other narrow Gaussian or by allowing the width of the sin- 



gle Gaussian to take a non-zero value [a ~ 200 eV). Since 
the parameters of these models are not well-constrained we 
prefer to present the results for the single Gaussian of zero 
width. (Also the other models imply continuum tempera- 
tures of ~ 2.2 keV which do not follow the trend of increas- 
ing continuum temperature with energy.) 

Although no specific calibration features at 7.8 keV have 
been reported we do mention that there exist some calibra- 
tion errors at high energies (i.e. above ~ 5 keV) due to the 
uncertainty in the optical axis, satellite wandering, and finite 
size of the Crab Nebula which is the key calibration source 
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Figure 14. As Figure but for Perseus. 



(see Gendreau & Yaqoob 1997). These warnings also apply 
to the line feature detected at 8.4 keV for Centaurus which, 
if of plasma origin, implies higher order K shell transitions 
from Fe (see Figure ^). 

Since the effects of calibration impact the various detec- 
tors differently, we have also measured the line properties for 
the SISO and SIS1 detectors separately. Generally we find ex- 
cellent agreement for the line properties between the SISO 
and SIS1. The line fluxes are usually slightly larger for the 
SISO, but the relative fluxes between the lines are generally 
consistent within the estimated 90 per cent errors for each 



detector. Although the SISO and SIS1 results are generally 
consistent, some features are easier to see in the SISO data; 
e.g. 8.4 keV feature of Centaurus. 

We have also verified the SIS results by simultaneously 
analyzing the SIS and GIS data. Only for the Fe K emis- 
sion is some tangible improvement in the constraints offered 
by adding the GIS data. The key improvement is in the 
continuum temperature which is consistent with the values 
determined from the SIS data but typically has about half 
the statistical error. However, because the SIS and GIS re- 
sponses differ we must let the line energies (as well as nor- 
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malizations) of both the SIS and GIS be free parameters in 
the fits (i.e. only the continuum temperature is tied between 
the detectors). The result is that adding the GIS does not 
appreciably improve the SIS constraints on the line fluxes. 

Although the GIS also does not detect the Fe XXVI 
(H) line near 6.9 keV in Centaurus and Perseus because it 
is blended with the Fe XXV complex, we do confirm the 
features near 8 keV in M87 and Centaurus with the GIS. 
The GIS does not resolve the 7.8 and 8.4 keV features in 
Centaurus, but the parameters of two gaussian features ob- 
tained at those energies with the GIS agrees with the SIS 
results. Similar to the SIS the 7.8 keV "feature" in Perseus 
is not obvious from visual examination of the GIS spectrum, 
though similar constraints on the feature are achieved with 
the GIS. The principal reasons why we present results for 
a 7.8 keV line in Perseus are (1) because such a feature is 
required in M87 and Centaurus, and (2) for comparison to 
previous work (see below). 

The results for the F XXV (He) line listed in Table | 
reflect a non-zero width on the Gaussian component. We 
obtained 90 per cent confidence a values of 48.20^p'gg 
eV for M87, 17.64±^;g eV for Centaurus, and 26.27^26.2? 
eV for Perseus; i.e. a is consistent with zero for Centau- 
rus and Perseus. When analyzed separately we found that 
a ~ 60(30) eV for the SISO (SIS1) data for the Fe line of 
M87. Similarly, a ~ 24(7) eV for the SISO (SIS1) data of 
Perseus. For Centaurus the best-fitting a ~ eV for each 
SIS. The significant a obtained for M87 may be due to finite 
intrinsic width of the Fe K complex, errors in the response 
matrix, or (most likely) a contribution to the emission from 
the (undetected) Fe XXVI (H) line. 

Our results represent the first detections of many of the 
line blends in these clusters. The most complete catalogue 
of X-r ay de tections was previously reported by Matsumoto 
et al. ( 1996] ) who used ASCA SIS data to constrain the Ka 
lines of H- and He- lik e Si an d S as well as Fe XXV (He) for 
M87. Fukazawa et al. ( 1994 ) presented equivalent widths for 
the Si XIV (H), S XV (He), S XVI (H), and Fe XXV (He) 
lines which are consistent with our resul ts (Ta ble Another 
notable paper is that of Molendi et al. (1998) who presented 
results for both the Ka and K/3 Fe XXV emission lin es in the 
center of Perseus using data from the SAX satellite ( Parmai 
et al. 19^) . Therefore we have presented the first constraints 
on (1) individual Ar and Ca lines in each of these systems, 

(2) the H-like Fe XXVI lines in Centaurus and Perseus, and 

(3) Fe K/3 emission in M87 and Centaurus. Our results for 
Mg XII (H) are also novel though our reported Mg flux is 
uncertain due to contamination from Fe L emission as noted 
above. However, it is reassuring that the energies of the fitted 
Mg lines lie exactly where Mg XII should be. 



3.2.2 Temperatures from line ratios and local continuum 

In Table ^| we list the temperatures implied by the H /He ra- 
tios of Ka lines of the same element for each element where 
90 per cent confidence limits on the fluxes are reported in 
Table Recall that for isothermal gas these ratios are very 
temperature sensitive (Figure |§]) and are independent of the 
elemental abundances. At the end of this section we shall 
also consider the temperatures suggested by the Fe K/3 lines. 
Since the temperatures listed in Table ^| are not obtained by 
fitting plasma models directly to the data they should be 



regarded as estimates which are sufficient for the qualitative 
discussion we present in this section. A quantitative assess- 
ment of the goodness-of-fit of different spectral models is 
given below in section 3.3. 



For M87 the plasma temperature indicated by the Si 
ratio is ~ 1.7 keV while the S ratio gives T ~ 2.1 keV. If the 
gas is isothermal then the temperature should be the same 
for the Si and S ratios. Instead the increase in temperature 
from the Si lines near ~ 2.0 keV to the S lines near ~ 2.5 
keV indicates that the emission measure consists of at least 
two components where the cooler component dominates the 
energy spectrum for E 2 keV while the hotter component 
dominates for higher energies. These results are also consis- 
tent with emission from a continuous range of temperatures 
where again the cooler components contribute most signifi- 
cantly to the Si ratios and the hotter components dominate 
the S ratio; note that for the integrated ASCA spectrum 
within the circular aperture a continuous range of temper- 
atures would be consistent with, e.g., the sum of different 
temperature components at different radii for a single-phase 
gas, or the sum of different phases in a multiphase cooling 
flow. 

The variation in local continuum temperature as a func- 
tion of energy (Table gives further support for multiple 
temperature components in the ASCA spectrum of M87; i.e. 
the continuum temperature varies from ~ 1.2 keV near Si 
to ~ 2.6 keV near Fe XXV. This variation in continuum 
temperature does not arise from an increasing proportion of 
bound-free emission since free-free emission dominates ex- 
cept near the strong Fe L lines near 1 keV (and for lower 
energies) . Although the local continuum temperature is con- 
sistent with the S ratio, it is somewhat lower than the Si 
ratio. This difference for Si could be due to an inaccurate 
continuum estimate arising from unresolved lines or from 
other systematics in the method we u se t o infer continuum 
temperature and line fluxes (section 3.2). However, for a 
multiphase plasma it is not necessary that the local contin- 
uum temperature agree with the temperature inferred from 
a line ratio assuming an isothermal gas. We also remark that 
the continuum temperature can be affected by non-thermal 
components whereas the line ratios are determined only by 
the hot plasma. 

Our results are not inconsistent with those obtained 



by Matsumoto et al. (1996) but a direct comparison of our 



results is not possible because they did not use a local con- 
tinuum when computing their line properties. These authors 
also examined the spectra within annuli of 2' width which 
may have been too small since the ASCA PSF is 3' in half- 
power diameter. Nevertheless, within our ~ 6' aperture Mat- 
sumoto et al. also found evidence for non-isothermal gas. 
They estimated that the continuum temperature increased 
from 2.1 to 2.3 keV within ~ 6' which is similar to the tem- 
perature difference we infer from the Si and S line ratios. 

The evidence for multiple temperature components is 
even stronger for Centaurus. The Ka line ratios imply tem- 
peratures of ~ 2.4 keV over energies ~ 2-3 keV, T ~ 3 keV 
near 4 keV, and T ~ 4 keV near 7 keV. The continuum tem- 
perature increases from T ~ 1.6 keV near 2 keV to T ~ 4 
keV near 7 keV. The good agreement between the temper- 
atures implied by the Fe line ratio and the continuum near 
7 keV indicates that a high-energy non-thermal component 
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Table 9. Temperatures Implied by ASCA Line Ratios 



Cluster 


Si 


S 


Ar 


Ca 


Fe 


Fe K/3/Ka 


M87 


i 7 +0.3 
i -' -0.3 


9 i+0.2 
z,1 -0.2 








Q + 5.8 

7 9 


Centaurus 


2 4+ ' 3 
Z - 4 -0.3 


9 9+O.3 
z - z -0.2 


24 +0.6 


1+0.7 
a,1 -l.l 


4 n+ - 4 




Perseus 


1 9+ - 4 

1 -0.2 




q Q+0:9 

°-°-0.9 


4 9+1.2 


5 9 +0 - 5 


2.8^-2 



Isothermal temperatures in units of keV corresponding to H /He flux ratios of the Ka lines listed in Table |8j also displayed are the Fe 
K/3/Kq ratios where Fe K/3 refers to the 7.8 keV feature and Fe Ka refers to the strong Fe XXV (He) line. For each line ratio we 
computed the (MEKAL) temperature appropriate for the resolution of the ASCA SIS as described in Figure [ll]. The error bars reflect 
90 per cent confidence limits. 



(e.g. AGN, discrete sources) does not dominate the contin- 
uum of Centaurus for energies ~ 5-9 keV. 

Hence, the variation of continuum temperature with en- 
ergy should reflect a related variation in the plasma tempera- 
ture, and when considering these temperatures and those in- 
ferred from the line ratios we conclude that the SIS spectrum 
of Centaurus clearly requires multiple temperature compo- 
nents for the hot gas within the ASCA aperture. (As with 
M87 the slight differences between continuum and line ratio 
temperatures are expected for a general spectrum with mul- 
tiple temperature component.) A variation in temperature 
between T ~ 2-4 keV is consistent with the two-temperature 
model s fitte d to the broad-band ASCA spectra by Fukazawa 
et al. < \l994j ). 

The Ka line ratios for Perseus give perhaps the 
strongest evidence for multitemperature hot gas of the three 
clusters in our sample. The ~ 2 keV temperature inferred 
from the Si ratio is significantly less than T ~ 6 keV im- 
plied by the Fe ratio. Similar to Centaurus the excellent 
agreement between the temperatures of the Fe ratio and the 
continuum near 7 keV argues against a strong contribution 
from a high-energy non-thermal component typical of, e.g. 
an AGN. In contrast to Centaurus the local continuum tem- 
peratures of Perseus are nearly constant with T ~ 5-6 keV 
over approximately 2-8 keV; the small glitch near S likely re- 
flects a miscalculation of the continuum due to our inability 
to compensate for the undetected He- like S XV line and the 
calibration e rror near 2.2 keV due to t he optical constants 
in the XRT ( |Gendreau fe Yaqoob 199^ . 

This near-isothermality of the continuum is not incon- 
sistent with the temperature variations implied by the Ka 
line ratios in Perseus. If the hottest temperature compo- 
nents near 6 keV dominate the continuum it is still possible 
for weaker cooler components to contribute to the line emis- 
sion at lower energies. For example, in Figure ^ we show that 
the HI Si XIV line is ~ 5 times stronger for T = 2 keV than 
for T = 6 keV for an isothermal gas; i.e. a cooler component 
can give a significant contribution to the line emission even 
if the continuum is dominated by hotter components. More- 
over, the cooler component contributes negligibly to the the 
Fe Ka ratio since the HI Fe XXVI line is over 100 times 
stronger for T = 6 keV than for T = 2 keV (Figure g). 
Thus, the Ka line ratios provide the key evidence for mul- 
titemperature gas in Perseus within the ASCA aperture. 

Finally, the ratios of K/3/Ka lines also probe the tem- 
perature of the gas (Figure |^) but since our estimates of the 
fluxes of the se lin es may be dominated by systematic errors 
(see section 3.2.1 ) we advise extra caution when interpret- 
ing Kp/Ka ratios. In particular, since the 2.2 keV features 
are pro bably dominated by systematic errors ( Sendreau & 
Yaqoot 1997), and the ratios of the fluxes of these features 



to those of the Si XIII (He) lines generally imply much larger 
temperatures than the other lines, we do not discuss these 
lines further. We also do not discuss the 1.7 keV features 
since (if of plasma origin) they are blends of lines of three 
different elements. 

The remaining candidate K/3 lines are the 7.8 keV fea- 
tures. As indicated in Table |^ the Fe K/3/Ka ratios are not 
very well constrained and in each case give temperatures 
consistent with the (H/He) Fe Ka ratio and/or the con- 
tinuum temperature near 7 keV; note that the ratio "Fe 
K/3/Ka" means the ratio of the flux of the 7.8 keV feature 
to that of the He-like Fe XXV line. For Centaurus the ratio 
of the flux of the 8.4 keV feature to that of the He-like Fe 
XXV line implies T ~ 9 keV for the best-fitting line ratio, 
but the 90 per cent confidence lower limit is T ~ 1.3 keV. 

Hence, these Fe K/3/Ka ratios do not require anoma- 
lously high temperatures at the 90 per cent confidence level 
and thus do not im plicat e resonance scattering as suggested 
by Molendi et al. ([1998) for Perseus from analysis of SAX 
data. Molendi et al. obtained an Fe K(3/Ka ratio of ~ 0.20 
within r — 6.4' which is about twice our value of ~ 0.10 
within r = 5.6'. However, the energy resolution of SAX is 
poorer than the ASCA SIS and thus to make a fair com- 
parison we need to add the H-like Fe XXVI flux to the Fe 
Ka emission; this reduces our Fe Kf3/Ka ratio to ~ 0.08 
corresponding to T ~ 1.7 keV. Although there may be sys- 
tematic errors in the ASCA data at these high energies as 
we discussed in section 



3.2.1 



(see Gendreau & Yaqoob 1997), 
the agreement we have have obtained between the contin- 
uum temperatures and the Fe Ka ratio for Perseus argues 
against a large systematic error. Therefore, we conclude that 
the Fe Kf3/Ka ratio in Perseus deduced from ASCA is con- 
sistent with optically th in pla sma and that the large ratio 
found by Molendi et al. (1998) with SAX must be due to an 
unrecognized systematic error. 

We mention that the Fe K/3/Ka results we have quoted 
assume a photospheric solar Ni/Fe ra tio. If instead we use 
the meteoritic solar ratios (see section 3.3) we find that the 
temperatures inferred for the Fe K/3/Ka ratios in Table ^| are 
reduced by a modest amount (« 25 per cent). To raise the in- 
ferred temperatures it is necessary to reduce the Ni/Fe ratio. 
Because the uncertainties are large we find that the temper- 
atures implied by the 7.8 keV/Fe Ka ratios are consistent 
with the local continuum temperatures even for Z-^i/Zp e ra- 
tios as small as ~ 1/5 (meteoritic solar). 



3.3 Broad-band spectral fitting 

The line ratios provide strong evidence for multiple tempera- 
ture components within the ASCA apertures in each cluster. 
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Table 10. Broad-Band ASCA Spectral Fits 



(keV) 
M87: 



(keV) 



EM C EM h 
(see notes) 



PL 

(see notes) 



Fc 



Si 



Ar 



Ca 



IT: (193.5/142/2.7e-3) 



2.16 



+0.04 



1.56"* 



1T+PL: (156.6/141/0.18) 



1 92+ - 07 



1.37 4 



2T: (130.3/140/0.71) 



1.06 



+0.23 
-0.42 



2.48 



+0.23 



0.36 



+0.22 



2T+PL: (125.8/139/0.78) 



0.85 



+0.31 



2.21 



+0.28 



0.22 



+0.20 



Constant £: (131.9/141/0.70) 



2.19^ 



2.54 



+0.24 



1.56 



+0.07 



1.25 



+0.16 
0.21 



1.29^ 



0.39 4 



0.20 



+0.15 



o.ss 1 



i.or 



0.78 



+0.13 



0.87 



+0.18 



0.75 



+0.13 



1 19+O.II n 97 +0.ie n r, +0.17 nal +0.31 

i ' 12 -0.10 U - 8 '-0.10 U - tli -0.17 U - 81 -0.30 

1 14+0- 14 94+ - 13 68+ - 21 1 12+ - 39 

• L " L *-a.l3 U - 94 -0.11 U - DO -0.19 i - iZ -0.37 



1 09+° 15 1 01+° 1 

1 - U3 _0.12 J - u — 0.1 

1 16+ ' 16 1 04+ ' 13 

1 - lu -0.14 A ' u *-0.12 



+0.12 n 00+O.21 n Qn +0.34 

i.ll U - b8 -0.20 U - 9U -0.33 

+0.13 n 79 +0.21 n qq+0.37 

-0.12 Ui ' -0.20 u - 3 -0.35 



104 +0.11 n 6 +0.H 
i-^-O.lO u - yD -o.ii 



+0.11 n 79 +0.21 n Q7 +0.36 
-0.11 u -' -0.20 u - 3 '-0.35 



Centaurus: 



IT: (489.1/301/3.9e-ll) 



3.23 



+0.05 



0.60 



+0.02 



1T+PL: (373.2/300/2. 6e-3) 



+0.10 



2.49 



2T: (295.0/300/0.57) 



0.39 



+0.04 



1.59 



+0.12 



4.80 



+0.39 



0.26 



+0.04 



2T+PL: (286.8/299/0.68) 

1 56 +0 17 4 09+ ' 64 21+ 006 
1 - OD -0.17 *' uy -0.56 u ' z± -0.07 

Constant £: (303.2/301/0.45) 

K + 0.06 5.Q2+ ' 24 1 RO+ 001 



3-46_ 06 
Perseus: 



0.62 1 



0.39 1 



0.33^ 



0.49"* 



0.22 



+0.10 
0.12 



1.28 1 



2.56"* 



1.24"* 



1.55"* 



1.26 



+0.08 



6.25 



+1.71 



1.75 



+2.09 
1.75 



1.76"* 



2.16 



1.31 



+0.27 



-0.09 
-0.09 



1.54 



1.39"* 



+0.20 



1.28 



+0.13 



1.73 



+0.24 



1.23 



+0.11 



1.46 



+0.22 



1.34 



+0.11 



0.68 1 



1.36" 1 



1.24? 



1.45 



l.W 



+0.30 



0.74^ 



1.50^ 



1.47"* 



1.62"* 



1.27" 1 



IT: (367.8/167/3.7e-17) 



4.62 



+0.07 



4.77 1 



1T+PL: (244.2/166/7.5e-5) 



+0.18 



3.21 



2T: (176.4/165/0.26) 



7 +0.i3 



,+0.29 



2 59+ ' 27 



1.32 1 



2T+PL: (168.8/164/0.38) 

' 1+0 ' 15 5.59+ ' 43 1 m+ - 19 



1.30"* 



1.13"* 



Constant £: (183.8/166/0.16) 

a+0.36 4 gl +0.08 



1.91 



9.03 1 



3.74 



+0.15 



3.05 



+0.42 



5.38 



+0.51 



2.00 



+ 1.09 
1.13 



0.62 



+0.03 



1.36 



0.70 1 



0.83 1 



+0.25 



0.76 



+0.04 



1.03^ 



1.44 



+0.23 



0.64 



+0.11 



0.77 



0.79^ 



+0.15 



0.00 



+0.07 



0.35 



+0.28 



0.51 



0.68 1 



0.37^ 



+0.26 



0.00 



+0.21 



0.59 



0.70^ 



0.89^ 



0.55^ 



+0.39 



Best-fitting and 90 per cent confidence limits on one interesting parameter (Ax 2 = 2.71) for models fitted to the summed SIS data 
between 1.6-9 keV. The models are defined similarly to those in Table Q with the following exceptions: (1) the values of \' 2 < dof, and the 
null hypothesis probability are listed in parentheses next to the model, (2) the emission measures (EM) are quoted in units of 
W~ 15 n e npV/4nD 2 similar to what is done in XSPEC, (3) PL indicates the normalization of the power-law model (photon index 1.7) 
at 1 keV i n uni ts of 10 — 2 photons cm -2 s — 1 keV -1 , and (4) the abundances are quoted in terms of the meteoritic solar values of 



Feldman (1992) and where "• ■ ■" indicates that the abundance of that element is tied to Fe in its meteoritic solar ratio. For models with 
multiple temperature components the abundances of each component are the same. Finally, the "Constant £" model represents a 
plasma emitting over a range of temperatures having equal emission measure over equal temperature intervals; for this model T c 
represents the temperature of the midpoint of the range while indicates the length of the temperature interval. 



To quantify more precisely the evidence for non-isothermal 
gas we must fit plasma models convolved with the detector 
response direc tly to the ASCA spectra. As we have shown 
in section 2.3 it is vital to consider as many line ratios as 



possible to distinguish between simple two-temperature and 
cooling flow models. Hence, in this section we examine mod- 



els fitted over 1.6-9 keV which includes the strong Si, S, Ar, 
Ca, and Fe K shell transitions. 

We do not analyze energies below 1.6 keV since the 
emission near 1 keV is dominated by Fe L emission, and in- 
trinsic absorbing material needs to be considered for lower 
energies (e.g. Fabian et al. 1994a). The latter point is es- 
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pecially relevant for our analysis because this allows us to 
assume (to a good approximation) that the same Galac- 
tic column density (Table ^) modifies each spectral compo- 
nent; i.e. this simplification limits the number of free pa- 
rameters for models that emit over a continuous range of 
temperatures. We account for absorption by our Galaxy us- 
ing the photo-electric absorption cr oss se ctions according to 
Balucihska-Church & McCammon (1992). 



Since it is our purpose to distinguish between rival spec- 
tral models we desire to optimize S/N by rebinning the 
ASCA sp ectra above the minimum 20 counts (see beginning 
of section 3.2) so that the energy resolution of the detectors 



are not greatly over-sampled. For the SIS data we settled 
on minimums of 50 counts for M87 and 100 counts for Cen- 
taurus and Perseus. Since the GIS resolution is poorer than 
for the SIS we chose larger minimums: 100 counts for M87 
and 150 counts for Centaurus and Perseus. The effects of 
using these larger PI bins are actually small, and the results 
we present below are consistent with the results obtained 
from using the standard 20-count binning employed in the 
previous section. 

We again focus on the SIS data to provide a consis- 
tent comparison with the results of the previous section, 
and because they provide the best constraints on the line 
properties. However, we discuss the additional constraints 
provide d by joint fitting of the SIS and GIS data below in 



section 3.3.3 



Unlike most previous studies we do not use th e pho - 
tospheric solar abundances of Anders & Grevesse (1989) 
but instead we use the meteoritic solar abundances of Fold- 



man (11993) which are more appropriate (Ishimaru & Ari- 
moto l f)97p . The principal difference is that the meteoritic 
Fe abundance (Fe/H of 3.24 x 10~ 5 by numb er) is 1/1.44 



times the photospheric value. Below in section 3.3.2 we dis- 
cuss how this choice affects our results in comparison to 
some previous studies. 



3.3.1 Constraints on f(T) 

In Table |l^ we list the results of fitting isothermal models 
(IT), two-temperature models (2T), and models with a con- 
stant differential emission measure, £(T) = constant, to the 
summed SIS data. The best-fittin g I T and 2T models are 
displayed at the bottom of Figures - [w]. The spike in the 
residuals of both the IT and 2T models of Perseus (Figure 
14|) coincides with the 2.2 keV feature discussed in section 
We find that this feature cannot be removed by mak- 



3.2.1 



ing simple alterations to £(T) which is consistent with the 
feature being dominated by the systematic error at 2.2 keV 
arising from the optical constants in the XRT ( Gendreau 
k. Yaqoob 1997). Hence, the results listed in Table |10| for 
Perseus exclude the energies 2.-2.8 keV where this feature 
dominates. 

From examination of the residuals in Figures p^l - 14] 
and the \ 2 nm l hypothesis probabilities (P) in Table |10 it 
is clear that IT models give poor fits (P <C 1) to the SIS 
spectra of each cluster. The most prominent residuals occur 
near the Si and Fe lines. We find that P decreases from 
M87 to Centaurus and then to Perseu s which agrees with 
our conclusion reached in section 3.2.2] using line ratios that 
an isothermal gas is most strongly excluded for Perseus and 



then Centaurus. The IT models for each cluster require a/Fe 
ratios different from solar to obtain a best fit. 

The 2T models provide superior fits that are formally 
acceptable (P > 0.1) for each cluster. Examination of Fig- 
ures [l2] - |l4| reveals no large residuals for these models 
though some excess emission at higher e nergie s may be in- 
dicated which we discuss below in section 3.3.S . Unlike M87 



the 2T models of Centaurus and Perseus are not improved 
significantly when allowing elements different from Fe to be 
varied separately: the 2T model having the abundances of 
the elements tied to Fe in their (meteoritic) solar ratios gives 
X 2 /dof/P values of 298.2/304/0.58 and 179.0/168/0.27 re- 
spectively for Centaurus and Perseus. (Note S is tied to Fe 
for Perseus since we exclude the energies 2-2.8 keV from the 
fits.) 

For comparison to the 2T model we examine the sim- 
plest model of a gas emitting over a range of temperatures 
with constant differential emission measure, £(T) = con- 
stant. This model provides fits that are virtually identical 
to the 2T models for M87 and Centaurus, and for Perseus 
the constant-^ model fits almost as well as the 2T model. 
Again the fits are not improved significantly for Centaurus 
and Perseus when allowing multiple abundances to be varied 
separately: we obtain 310.1/305/0.41 and 195.4/169/0.08 
respectively for Centaurus and Perseus for the constant-^ 
model where only the total metallicity is varied. We em- 
phasize that these constant-^ models have qualitatively dif- 
ferent emission measure distributions from the 2T models. 
That is, for lower and upper temperatures denoted by Tl ow 
and Tiii, the constant-^ models indicate that the hot gas 
emits at constant emission measure over the following tem- 
perature ranges: Ti ow « 0.92 keV and Tm ~ 3.46 keV for 
M87; Ti ow » 0.92 keV and T hi « 5.9 keV for Centaurus; 
Ti ow w 0.43 keV and T hi ta 9.4 keV for Perseus. 

The strong similarity between the fits of the 2T and 
constant-^ models highlights the difficulty in constraining 
the shape of a general £(T) with ASCA data. Although we 
have not presented re sults for specific multiphase cooling 
flow models (sections 2.2.3 and |2.3| ) because of computa- 
tional expense, the similarity between the fits of the 2T and 
constant-^ models is sufficient to demonstrate that multi- 
phase cooling flows also cannot be distinguished from the 
2T models; i.e. multiphase cooling flows have intermedi- 
ate between the 2T and constant-^ models. Therefore, the 
fits to the SIS data over 1.6-9 keV for each cluster clearly 
rule out isothermal gas but cannot determine whether the 
gas emits at only two temperatures or over a continuous 
range of temperatures as expected for a multiphase cooling 
flow. 



3.3.2 Abundances 

The Fe abundances obtained from isothermal models are 
very similar to those obtained from the multitemperature 
models, whereas several of the a/Fe ratios differ from solar 
for each IT model but not for the multitemperature models 
of Centaurus and Perseus. (A possible exception being the 
sub-solar Ar/Fe ratio for the constant-^ model of Perseus.) 
The fact that the Si/Fe ratios exceed solar for IT models in 
contrast to the 2T and constant-^ models is another demon- 
stration of biased measurements of abundances which occur 
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from fitting a multitemperature spectrum with an isother- 
mal model (e.g. see Appendix of Buote 1999b). 

If instead the photospheric solar abundances are used 
as done in most previous studies, then a/Fe ratios in ex- 
cess of solar are required to obtain fits of similar quality to 
the (meteoritic) multitemperature models of Centaurus and 
Perseus; e.g. using photospheric solar abundances for the 2T 
model of Centaurus we obtain Z-p c — 0.83Zq Zsi = 1.2Zq 
Z s = 1-IZq Zkt = IAZq Z Cs , = 1.4Z . This Si/Fe ra- 
tio of ~ 1.5Zq (photospheric) is similar to the value of 2 
solar (photospheric) obtained by Fukazawa et al. ( 1998 ) us- 
ing isothermal models within an annulus of 0.1/ig 



0.4/i 



Mpc which excludes the region we have analyzed (Table | 
When using the photospheric solar abundances for the IT 
model of Centaurus we obtain a Si/Fe ratio of 2 solar in 
good agreement with Fukazawa et al.. Hence, the Si/Fe en- 
hancements obtained by Fukazawa et al. are due in part to 
their using photospheric solar abundances, and if the hot gas 
within their apertures deviates significantly from isothermal 
then the assumption of isothermal gas contributes as well. 

When accounting for the differences between meteoritic 
and photospheric solar abundances we find that the Fe abun- 
dances we have obtained agree well with previous ASCA 
studies; e.g. for M87 (Matsumoto et al. 1996; Hwang et al. 
1997), for Centaurus (Fukazawa et al. 1994; Fabian et al. 
1994a), and for Perseus (Fabian et al. 1994a). We mention 
that for the multitemperature models we do not require the 
abundances to be different for different temperature compo- 
nents. That is, our models do not require abundance gradi- 
ents within our apertures but they do not exclude them ei- 
ther. The existence of abundance gradients within our aper- 
tures does not compromise the abundances inferred by our 
simple models because, as we have shown previously for el- 
lipticals and groups (see section 4.1 in Buote 1999b and 
section 5.3 of Buote 1999a), these models provide accurate 
measurements of the average Fe abundances obtained from 
the integrated spectrum within a large spatial region. 

3.3.3 Evidence for high energy excesses? 

Examination of the residuals from the 2T models (Figures 
|l2] - |l^) suggests excess emission for energies above ~ 7 
keV for Centaurus and perhaps the other clusters as well. 
Using simple models we have investigated the contribution 
of a non-thermal component which manifests itself preferen- 
tially at higher energies in the ASCA bandpass. We find that 
the results for a canonical AGN component (i.e. power law 
with photon index 1.7 - Mushotzky et al. 1993) are nearly 
identical to those of a bremsstrahlung component with tem- 
perature greater than ~ 10 keV. We shall focus on the results 
using the AGN model. 

In Table [l(] we give the results of adding the simple 
AGN power-law model (PL) to the IT and 2T models. The 
IT models are improved significantly when the PL compo- 
nent is added, but the improvement in each case is not nearly 
as great as that observed for the 2T models over the IT mod- 
els. The PL component only improves the 2T models slightly. 
A negligible improvement in the multitemperature fits is ex- 
pected because of the good agreement between the Fe Kq 
ratios and th e local continuum temperatures discussed in 



eluded in the fits then, because of its better sensitivity to 
energies above ~ 8 keV, the need for a power-law component 
will be made clearer. In fact, we find just the opposite. When 
the GIS data are included the relative improvement of the 
2T+PL model over the 2T models is Jess for Centaurus and 
Perseus although for M87 the results are similar to those 
of the SIS data alone. We find that for the simultaneous 
SIS+GIS fits for Centaurus and Perseus the SIS residuals 
exceed the 2T model while the GIS residuals lie below the 
best-fitting 2T model. Differences between the SIS and GIS 
at these energies are not unexpect ed be cause of the cali- 
bration issues mentioned in section 3.2.1 (see Gendreau & 
Yaqoob 1997). At any rate, the weak evidence for the power- 
law components obtained from the 2T models (which also 
applies t o cons tant-^ models) is unclear because Gendreau & 
Yaqoob (1997) note that existing calibration errors can give 
rise to a hard tail in the data. Hence, our spectral analysis 
does not provide convincing evidence in support of excess 
hard emission over that produced by the hot plasma compo- 
nents. The lack of excess hard emission in M87 is sup ported 
by recent RXTE observations ( Reynolds et al. 1999[ ) . 

We mention as a caveat that if lower energies are in- 
cluded in the fits (e.g. down to ~ 0.6 keV) then the temper- 
atures of the 2T models are lowered (due to the influence 
of cooler temperature components) and there is more of a 
need for an additional higher energy component which can 
be modeled with a (now intrinsically absorbed) power law. 
However, disentangling the uncertain shape of £(T) and the 
relative contribution of a power-law component from issues 
of excess absorption and possible uncertainties in the Fe L 
lines is beyond the scope of this paper; i.e. we conclude that 
the multitemperature models fitted over 1.6-9 keV, which 
do not require us to accurately model excess absorption or 
the Fe L lines, also do not require high energy excesses. 



4 PROSPECTS FOR FUTURE MISSIONS 

ASCA has provided a tantalizing glimpse of the temperature 
structure from emission line ratios in a few ellipticals and 
clusters. We are fortunate that within the next few years 
satellites will provide superior X-ray data capable of deter- 
mining £(T) with unprecedented precision. Our intention in 
this section is to give some indication of how well the multi- 
phase structure of ellipticals and clusters can be probed by 
these new missions. 

We continue with our prototype examples of NGC 4472 
and Centaurus to examine how well these new missions can 
distinguish a co oling flow from a two-temperature plasma 
(see Section 2.3). That is, we simulate observations of cool- 
ing flow models of these clusters (Table ^) and then fit two- 
temperature models to the simulated data. We then compare 
ratios of the best-fitting two-temperature models to the sim- 
ulated data. 



4.1 Chandra and XMM 



Our probe into the future begins with the Chandra X-Ray 
Observatory]^] which is scheduled to launch in 1999 and the 



section 3.2.2 



It might be expected that when the GIS data are in- 



tt 
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Figure 15. Simulated 100 ks observations of NGC 4472 and Centaurus with the Chandra ACIS-I following the cooling flow models in 
Table W The best-fitting two-temperature models are also shown. 
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channel energy (keV) channel energy (keV) 

Figure 16. Same as Figure ^ except we plot the ratio of the best-fitting two-temperature model to the simulated data arising from 
the cooling flow model. Note that we increased the energy bin sizes with respect to Figure ha for aesthetic reasons. 



X-Ray Multi-Mirror Mission (XMM^\ which is scheduled 
to launch in January 2000. The ccds on board these satellites 
have similar energy resolution, and thus we shall focus on 
simulations of Chandra observations and note any important 
differences expected from XMM. 

There are several instruments on board Chandra, but 
we focus our attention on the AXAF CCD Imaging Spec- 
trometer (ACIS) since it is most appropriate for high reso- 
lution spectral analysis of extended sources. The energy res- 
olution (and bandwidth) of the four front illuminated CCDs 
employed in the ACIS imaging array (ACIS-I) is similar to 
the pre-fiight ASCA SIS: e.g. AE ~ 60 eV (FWHM) at 
E ~ 1 keV and AE ~ 140 eV (FWHM) at E ~ 6.5 keV. 
However, the spatial resolution of the ACIS-I is vastly supe- 
rior: 90 per cent encircled energy fraction within 1" of a 1.5 
keV monochromatic point source. 

With XSPEC we simulated 100 ks ACIS-I observations 
of NGC 4472 and Centaurus using the cooling flow mod- 
els listed in Table §; we used the default RMF and ARF 
files supplied for simulating AOl observations. We included 
background in these simulations using models for the diffuse 



X-ray background and instrumental background following 
Figure 5.14 of the AXAF Proposers' Guide Rev. l.oQ. The 
background-subtracted spectrum of NGC 4472 is shown in 
Figure [Hj. Also shown is the best-fitting 2T+BREM model 
(variable A% for each plasma component, abundances are 
tied together). In Figure |l6| we plot the ratio of the best- 
fitting 2T+BREM model to the simulated data. 

The simulated ACIS-I data for NGC 4472 are of sub- 
stantially higher quality than the real ASCA data (see Buote 
1999a). Line blends of O VIII (H), Mg XI (He) and XII (H), 
Si XIII (He) and XIV (H), S XV (He) and XVI (He), and 
Ar XVII (He) are readily apparent. Deviations of the model 
and simulated data are apparent for the Ka lines of Si and 
S, but the discrepancies are only significant at ~ 90 per cent 
level. However, the 2T+BREM model is clea rly r ejected by 
lines at lower energies. As expected (section 2.3) the most 
important residuals occur for the Fe L lines from ~ 0.7 — 0.9 
keV and the O VIII (HI) line. 

We note that XMM is able to reject the 2T+BREM 
model with higher confidence than Chandra. Firstly, since 
XMM has considerably larger effective area the constraints 



M http://astro.estec.esa.nl/XMM/ 
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Figure 17. Same as Figure ^ except the observations are simulated for the ASTRO-E XRS. 
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Figure 18. Same as Figure ^ except we plot the ratio of the best-fit two-temperature model to the simulated data arising from the 
cooling flow model. Note that we increased the energy bin sizes with respect to Figure |l^ for aesthetic reasons. 



on the Ka lines of Si and S are more precise; i.e. it is 
not so necessary to appeal to the Fe L lines to reject 
the two-temperature model. Secondly, the larger bandwidth 
(0.1-12 keV) of XMM makes it more difficult for the two- 
temperature model to mimic the cooling flow spectrum. 

The simulated 100 ks observation of Centaurus and the 
best-fitting 2T model are also shown in Figure |l5|; similarly 
the ratio of the best-fitting 2T model to the simulated data is 
plotted in Figure All of the line blends detected with the 
ASCA data (Figure ^) are also clearly detected by Chan- 
dra. Similar to NGC 4472, the line fluxes are better con- 
strained than with ASCA. The H/He Ka ratios of Si and S 
are both marginally inconsistent with the two temperature 
models at the ~ 2a level as is apparent from the excesses in 
Figure [h| Again the most substantial deviations of the two- 
temperature model occur for energies ^ 1 keV where the Fe 
L and oxygen lines dominate the spectrum and easily reject 
the two-temperature fit. Nevertheless, it is significant that 
the H/He Ka ratios of Si, S, and Fe can independently dif- 
ferentiate the two-temperature and cooling flow models for 
Centaurus so that issues like the accuracy of the Fe L tran- 
sitions in the plasma codes and the amount of absorption 
are not an issue. 



4.2 ASTRO-E 

After Chandra and XMM a significant leap forward in en- 
ergy resolution will be provided by the X-Ray Spectrom- 
eter (XRS) on board the ASTRO-E satellite^ scheduled 
for launch in March 2000. Unlike previous X-ray detectors, 
the XRS uses an array of X-ray micro-calorimeters which 
will provide an unprecedented combination of high energy 
resolution (AE ~ 10 eV FWHM) and large collecting area. 
Although its spatial resolution is limited and its expected 
lifetime is only ~ 2 years, the XRS will have the ability to 
measure individual emission lines of extended sources such 
as ellipticals, groups, and clusters like no previous X-ray 
mission. 

In analogy to the Chandra cases, we simulated 100 
ks observations with the XRS for the cooling flow models 
of NGC 4472 and Centaurus. We used the default RMF, 
ARF, and background files provided by the ASTRO-E GOF 

web sitejnjl We show in Figure [l?] the simulated XRS spec- 
tra and best-fitting two-temperature models of NGC 4472 
(2T+BREM) and Centaurus (2T). The ratios of the two- 



11 http: / /lheawww. gsfc . nasa.gov /docs/ xray / astroe / 
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temperature models to the simulated data are displayed in 
Figure |l| 

With the vastly improved resolution many blends of 
lines in the Chandra and XMM (and ASCA) data separate. 
The He- like Ka transitions of Si and S become separate He4 
lines and blends of He4-5 and satellites; similar transitions 
of Mg are also resolved. Many new transitions are resolved 
including the K/3 lines of Si XIII (He3) and Si XIV (H2). 
The Fe XXV complex for Centaurus, though still unresolved, 
shows significant structure indicating separate contributions 
from the He4, He5-6, and satellites. Large residuals in the 
two-temperature fits at the energies of these K shell lines 
shows that finally at the resolution of the XRS the K shell 
lines can clearly differentiate between two-temperature and 
cooling flow models for NGC 4472 in addition to Centaurus. 
The novel constraints on ratios of K/3/Ka lines of the same 
clement can be used to verify fundamental assumptions of 
the coronal model; e.g. the assumptions of ionization equi- 
librium and an optically thin plasma. 

Moreover, the oxygen lines and many of the key Fe L 
blends are easily measured with the XRS. The O VII (He) 
and OVIII (HI) lines for NGC 4472 clearly distinguish the 
simulated cooling flow from the two-temperature fit. (Note 
that allowing the relative abundances of, e.g., O, Ne, Mg, Si, 
S, Ar, and Ca to be free does not improve the fit substan- 
tially.) Of more importance is the Fe XVII blend at ~ 0.7 
keV (see Table ^|) and other Fe L lines near 1 keV. It is clear 
that data of this quality will provide exciting new insights 
into the temperature structure of the hot plasma in these 
systems. 



5 DISCUSSION 

5.1 Resonance Scattering 

It is usually assumed that the hot gas in ellipticals and clus- 
ters is optically thi n to X-rays. However, Gil'fanov, Syun- 
yaev, & Churazov (1987) have argued that resonance scat- 
tering of photons produced by strong transitions to the 
ground state of an abundant ion may be important in these 
systems. Gil'fanov et al. estimated significant optical depths 
r ~ 4.0 and 2.7 to the Fe XXV (He4) line in the centers of 
Perseus and M87 respectively. They predicted that steeply 
rising abundance profiles with increasing radius should be 
observed if the intrinsic abundances are constant with ra- 
dius. 

In fact, observations of clusters almost always show flat 
abundance profiles or, particularly for cooling flow clusters, 
abundance profiles that increase towards the center (e.g. 
Fukazawa et al. 1994). It would seem unlikely that in most 
cases the nearly flat abundance profiles of observed clusters 
are actually due to large intrinsic negative abundance gradi- 
ents which almost exactly cancel the predicted effect due to 
resonance scattering. As noted by Gil'fanov et al., turbulent 
motions in the gas may significantly reduce the importance 
of resonance scattering. (Turbulent velocities of ~ 300 km 



As shown in section 3.2.2 the Fe K/3/Ka ratio measured 



from the ASCA data of Perseus is consistent with optically 
thin plasma and does not suggest resonance sca tterin g in 
contrast to the result obtained by Molendi et al. (1998) us- 



s should be sufficient to reduce the optical depth of the 
Fe XXV He4 line below unity in a massive cluster.) For- 
tunately, as shown by Gil'fanov et al., resonance scattering 
modifies the observed line shapes and thus its importance 
may be directly measured by future missions with sufficient 
energy resolution (e.g. Constellation-X) . 



ing SAX data. Although resolving the discrepancy between 
the ASCA and SAX results is beyond the scope of our pa- 
per, we believe the consistent results we have obtained for 
M87, Centaurus, and Perseus using ASCA data argue that 
resonance scattering has a negligible impact on the Fe abun- 
dances inferred from current X-ray data. 

5.2 Multiphase Strengths of Clusters and Qq 

Since structure formation in the Universe is believed to pro- 
ceed by hierarchical clustering, galaxy clusters should show 
evidence of merging depending on their dynamical "age". 
The relative fraction of dynamically young clusters in a com- 
plete cluster sample should be sensitive to the cosmological 
density pa ramet er Qo as first suggested by Richstone, Loeb, 
& Turner ( 1992 ) . Recent analysis of X-ray cluster morpholo- 
gies favors a unive rse with low matter density fio ~ 0.35 
( Buote fe Xu 1997j ). 

Cluster morphologies, however, are only one indicator 
of the dynamical states of clusters. The dynamical state in- 
dicated by cluster morphology is the same as that indicated 
by the quantity, AM/ M, the fractional amount of mass a 
cluster accreted during the previous crossing time, but the 
connection between cluster morphology and AM/ M (and 
thus Qo) is rigorously valid only for clusters that are not very 
close to a virialized state (see section 2.1 of Buote 1998). 

In contrast, the temperature structure of clusters can 
provide a sensitive measure of the dynamical states of nearly 
relaxed clusters. For illustrative purposes, let us consider a 
massive cluster (T ~ 8 keV) which has recently formed. At 
this time there will be little or no contribution to the emis- 
s ion measure (£(T)dT) fro m a cooling flow since arguments 
( |McGlynn fe Fabian 1984] ) and ROSAT data (Buote & Tsai 
1996; see also Jones & Forman 1999) show that the pres- 
ence of cooling flows anti-correlates with substantial sub- 
structure. Within the central few hundred kpc there will be 
significant temperature fluctuations arising from incomplete 
relaxation (e.g. as in A754 and A2256 - Henry & Briel 1995 
and Briel & Henry 1995). When the cluster has relaxed to a 
state that is close to the onset of a cooling flow, the emission 
measure distribution will be peaked near T ~ 8 keV along 
with residual temperatur e fluc tuations. In terms of the mul- 
tiphase strength (Section 2.2.1), such a cluster about to start 
a cooling flow will have typically cr^ 0.1. 

As the cluster relaxes a cooling flow will gradually de- 
velop. The emission measure distribution will widen until it 
resembles those in Figure ^ with <tj ~ 0.7 - 0.8. (Within 
the central few hundred kpc the cooling flow can contribute 
up to ~ 80 - 90 per cent of the X-ray emission; e.g. Fabian 
1994.) Notice that the dynamic range in between a cool- 
ing flow cluster (~ 0.7) and an incipient cooling flow 0.1) 
is considerably larger than the difference between an incip- 
ient cooling flow and an isothermal gas (crj = 0); i.e. the 
temperature structure, expressed in terms of er^, is a sen- 
sitive discriminant of the "age" of nearly relaxed clusters 
unlike cluster morphologies. 

However, clusters which have recently experienced ma- 
jor mergers can have a wide distribution of temperature 



© 0000 RAS, MNRAS 000, 000-000 



28 D. A. Buote et al. 



components contributing to the emission measure, which in 
some cases may give crj comparable to cooling-flow clus- 
ters. It will be necessary, as a result, to consider the joint 
(Tf and morphology (e.g. log 10 P m /Po - Buote & Tsai 1995) 
distributions to remove this potential degeneracy of clus- 
ter ages with large o^. Hence, with future X-ray satellites 
capable of detailed mapping of the emission measure dis- 
tributions of clusters (e.g. ASTRO-E, Constellation-X) , the 
multiphase strength can be measured for many clusters 
and be compared to those produced in N-body/ Hydrody- 
namical simulations (and combined with information from 
cluster morphologies) in order to constrain fio- 

It is possible with currently available data from ASCA 
to estimate cr^ of many clusters without detailed mapping of 
£ (T) for each cluster. This is possible if one assumes that the 
true £(T) of clusters is accurately represented by a two com- 
ponent model consisting of a constant pressure cooling flow 
and an isothermal ambient gas component (e .g. co oling flow 
models in Table ^J). As described in section 2.2.1, the mul- 
tiphase strength of such models is just /oj , where er^ is the 
multiphase strength of the cooling flow component (Figure 
^|) and / is the fractional contribution of the cooling flow 
to the total emission measure. Since these two-component 
models describe the available data of clusters quite well, this 
procedure for estimating cr^ should be a good approximation 
in most cases. 



6 CONCLUSIONS 

We have studied the ability of X-ray emission line ratios to 
probe the temperature structure of elliptical galaxies and 
galaxy clusters using theoretical models, the best available 
ASCA data, and simulated observations with Chandra and 
ASTRO-E. We have emphasized multiphase spectra consist- 
ing of either two temperatures or cooling flow models and 
examined the best means of distinguishing between them. 
Ratios of lines of K shell transitions have been emphasized 
over the Fe L lines because of possible remaining inaccura- 
cies in the available plasma codes. 

From analysis of isothermal MEKAL models computed 
with energy resolution 2 eV (FWHM) appropriate for the 
proposed Constellation-X mission, we found that the most 
useful ratios of K shell transitions for probing the emis- 
sion measure distributions of ellipticals and clusters are the 
Hl/He4 (Kq) and satellite/He4 ratios of the same element 
(notation described in Table §). The H2/H1 and He3/He4 
ratios (K/3/Ka) are much less temperature sensitive for tem- 
peratures where the transitions are strong. Ratios of similar 
transitions of different elements (e.g., He4 of S/Fe) are very 
temperature sensitive but must be used with caution for in- 
terpreting temperatures structure since they depend on the 
relative abundances of the elements. Finally, we also discuss 
the very strong temperature dependence of a few ratios of 
Fe L lines between ~ 0.7 - 1.4 keV which are typically the 
strongest lines and most temperature sensitive ratios of all. 

We considered a few simple multiphase models to evalu- 
ate how well the various line ratios can constrain the general 
differential emission measure £(T). To facilitate quantita- 
tive comparison of ellipticals and clusters having different 
£(T) we introduced the "multiphase strength" parameter, 
<T(, which measures the fractional width of f(T): ctj = 



for an isothermal gas and crj = 1 for a constant emission 
measure distributed over a large temperature range. Of the 
simple multiphase models examined (two temperature, con- 
stant emission measure, gaussian emission measure, cooling 
flows) it is found that the two-temperature models have the 
most flexibility; i.e. for a given value of (T) (emission mea- 
sure weighted temperature) the two-temperature models can 
generally provide the best approximation to a given line ra- 
tio. 

The flexibility of the two-temperature model is illus- 
trated by its ability to match very closely a model consisting 
of a cooling flow component and a single temperature com- 
ponent corresponding to ambient gas (Figure p"o| ). Cooling 
flows have nearly maximal multiphase strength (er^ ~ 0.7), 
but their Ka line ratios (aside from O) can be approximated 
by two temperature models within ~ 20 per cent for ellip- 
tical galaxies like NGC 4472 and within ~ 50 per cent for 
clusters like Centaurus. The best means to distinguish be- 
tween two-temperature models and cooling flows is to use 
oxygen lines or Fe L lines whose ratios can differ by a factor 
of ~ 10 between the models. 

We have re-analyzed the ASCA spectra of three of the 
brightest galaxy clusters to assess the evidence for multi- 
phase gas in their cores: M87 (Virgo), the Centaurus clus- 
ter, and the Perseus cluster. Ka emission line blends of Si, 
S, Ar, Ca, and Fe are detected in each system as is signif- 
icant Fe K/3 emission. The Fe K/3/Ka ratios are consistent 
with optically thin plasma models and do not suggest res- 
onance scattering in these systems. Consideration of both 
the ratios of H-like to He-like Ka lines and the continuum 
temperatures clearly rules out isothermal gas in each case. 

To obtain more detailed constraints we fitted plasma 
models over 1.6-9 keV where the emission is dominated by 
these K shell lines and by continuum. In each case the ASCA 
spectra cannot determine whether the gas emits at only two 
temperatures or over a continuous range of temperatures as 
expected for a cooling flow. The metal abundances are near 
solar for all of the multiphase models. For Centaurus and 
Perseus the Si/Fe ratios are consistent with the meteoritic 
solar values. We also find that our multitemperature models 
fitted over 1.6-9 keV, which do not require us to accurately 
model excess absorption or the Fe L lines, also do not require 
any excess emission at higher energies from, e.g, an AGN 
component. 

We examined the ability of future X-ray missions to 
constrain £(T) in ellipticals and clusters. By simulating ob- 
servations with the Chandra ACIS-I and the ASTRO-E XRS 
of cooling flow models of NGC 4472 and Centaurus we found 
that these future instruments will have the ability to distin- 
guish between two-temperature models and cooling flows. 
However, the lower resolution Chandra ACIS-I (and also the 
XMM ccds) must rely on constraints from O and Fe L lines 
for definitive results whereas the ASTRO-E XRS only re- 
quires the Ka lines above ~ 2 keV (e.g. Si and S). 

We have discussed the evidence for resonance scattering 
in light of our results for the K/3/Ka ratios. We also discuss 
the multiphase strength, crj, as an indicator of the dynamical 
states of clusters and its potential use as a probe of Qo in 
conjunction with galaxy cluster morphologies. 

The difficulty in distinguishing between different multi- 
phase models, specifically a two-temperature plasma from a 
cooling flow, has been emphasized in this paper. Presently, 
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there are only a few cases where line ratios have been mea- 
sured with sufficient precision to merely indicate that the 
gas is not isothermal. It is imperative that the temperature 
structures of ellipticals and clusters be considered when in- 
terpreting properties (e.g. abundances) derived from spec- 
tral fits in these systems. 
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